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We investigate the LHC discovery potential of _R-parity violating supersymmetric models with a 
right-handed selectron or smuon as the lightest supersymmetric particle (LSP). These LSPs arise 
naturally in _R-parity violating minimal supergravity models. We classify the hadron collider signa- 
tures and perform for the Hrst time within these models a detailed signal over background analysis. 
We develop an inclusive three-lepton search and give prospects for a discovery at a center-of-mass 
energy of ^/s = 7 TeV as well as ^/s — 14 TeV. There are extensive parameter regions which the LHC 
can already test with y^s = 7 TeV and an integrated luminosity of 1 fb~^. We also propose a method 
for the mass reconstruction of the supersymmetric particles within our models at — 14 TeV. 



I. INTRODUCTION 

Since 2010, the Large Hadron Collider (LHC) is col- 
lecting data at a center of mass energy of -^s = 7 TeV 
and first searches for physics beyond the Standard 
Model (SM) have been published [IHll. Even with 
only an integrated luminosity of 35 pb~^, the LHC 
has already tested supersymmetric models [11, be- 
yond the Tevatron searches [3 ■ Furthermore, it 
is expected that the LHC will collect 1 fb~^ of data 
until the end of 2011. 

One of the most promising LHC signatures for su- 
persymmetry (SUSY) are multi-lepton final states [iB - 
hj. On the one hand, electrons and muons are easy 
to identify in the detectors. On the other hand, the 
SM background for multi-lepton final states is low. In 
this publication, we focus on such signatures. 

We consider the supersymmetric extension of the 
SM with minimal particle content (SSM) [H, EM- 
Without further assumptions, the proton usually has 
a short lifetime in this model [18M20|, in contradic- 



tion with experimental observations |21| . The pro- 
ton decays, because renormalizable lepton and baryon 
number violating interactions are jointly present. One 
therefore must impose an additional discrete symme- 
try. The most common choice for this discrete symme- 
try is i?-parity, or equivalently at low-energy: proton- 
hexality (Pg)- Either suppresses al l lep ton- and baryon 
number violating interactions (22H24j . The SSM with 
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i?-parity is usually denoted the minimal supersymmet- 
ric SM (MSSM). 

We consider here a different discrete symmetry, 
baryon-triality (B3) (22| - [25j . which suppresses only the 
baryon number violating terms, but allows for lepton 
number violating interactions. The B3 SSM has the 
advantage that neutrino masses are generated nat- 
urally [26|-|29[ without the need to introduce a new 
(see-saw) energy scale j30l - l32l |. The lepton number 
violating interactions can be adjusted, such that the 
observed neutrino masses and mixing angles can be ex- 
plained [H,!!^]- Note that both Pk and B.^ are discrete 
gauge anomaly free symmetries j2a - l24l . l3a | . 

In the B3 SSM, the lightest supersymmetric parti- 
cle (LSP) will decay via the lepton number violating 
interactions and is thus not bounded by cosmological 
observations to be the lightest neutralino, Xi Un- 
like in the MSSM, the Xi is not a valid dark matter 
(DM) candidate. However, several possible DM can- 
didates are easily found in simple extensions of the 
B3 SSM; for example, the axino |38l - l4l| . the gravitino 
[^Hl] or the lightest [/-parity particle [4i.l45|. 

We consider in this paper the B3 SSM with a right- 
handed scalar electron (selectron, cr) or scalar muon 
(smuon, flfi) as the LSP. These LSP candidates natu- 
rally arise in the B3 minimal supergravity (mSUGRA) 
model [4^, on which we focus in the following. Here, 
large lepton number violating interactions at the grand 
unification (GUT) scale drive the selectron or smuon 
mass towards small values at the electroweak scale via 
the renormalization group equations (RGEs) [l^l ■ We 
describe this effect and the selectron and smuon LSP 
parameter space in the next section in more detail. 
Further LSP candidates within B3 mSUGRA (beside 
the Xi) are the lightest stau, fi 0, IH, H^, and the 
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sneutrino, i^e,fj.,T [46l l49j. depending on the dominant 
lepton number violating operator j47| . 

If SUSY exists, the pair production of strongly inter- 
acting SUSY particles (sparticles), like scalar quarks 
(squarks), is usually the main source for SUSY events 
at hadron colliders like the LHC [s^l- Furthermore, 
squarks, are much heavier than the Xi in most su- 



persymmetric models [SlJ. Assuming that we have a 
right-handed selectron or smuon, in, as the LSP, a 
natural cascade process at the LHC is 



(1) 



where the squarks decay into a quark, q, and the Xi- 
The Xi decays into the in LSP and an oppositely 
charged lepton, i, of the same flavor. 

The Ir LSP can then decay via the lepton number 
violating interactions, for example 



as a function of the dominant i?-parity violating inter- 
action. Based on this, we develop in Sec. |IV] a set of 
cuts for an inclusive three-lepton search at the LHC 
and give prospects for a discovery at ^/s = 7 TeV 
as well as at ^/s ~ 14 TeV. In Sec. |V] we propose a 
method for the reconstruction of the supersymmetric 
particle masses. We conclude in Sec. IVIl 

Appendix [previews the mass spectrum and branch- 
ing ratios of our benchmark models and Appendix |B] 
shows the cutflow for our y/s — 14 TeV analysis. We 
give in Appendix [C] the relevant equations for the kine- 
matic endpoints for the mass reconstruction of Sec. fVl 
and calculate in Appendix ID] some missing 3-body de- 
cays of sleptons. 



II. THE SELECTRON AND SMUON AS THE 
LSP IN R-PARITY VIOLATING MSUGRA 



The B3 mSUGRA Model 



i.e. into another charged lepton £' and a neutrino ly. 
As we argue in the following, this is the case for large 
regions of the B3 SSM parameter space. We thus ob- 
tain from Eqs. ([T|) and ([2]) an event with four charged 
leptons in the final state. Taking into account that 
some leptons might not be well identified, we design 
in this paper an inclusive three-lepton search for in- 
LSP scenarios. Although we concentrate on the B3 
mSUGRA model, our results apply also to more gen- 
eral models as long as Eqs. ((T)) and © hold. We will 
show that because of the high lepton multiplicity in B3 
models, the discovery reach at the LHC with ^/s = 7 
TeV exceeds searches in the i?-parity conserving case 
[53 |. We also give prospects for a discovery at ^/s = 14 
TeV and propose a method for the reconstruction of 
sparticle masses within our model. 

The phenomenology of slepton LSPs has mainly 
been investigated for the case of a stau LSP. See for ex- 
ample Refs. [ll,|3l,|4i,lil[H-|6l. Recently, Ref. ^ 
proposed a tri-lepton search for stau LSP scenarios, 
which is similar to our analysis, although the stau 
in Ref. [ll] decays via 4-b ody decays. LEP II has 
searched for slepton LSPs [gj, ■ No signals were 
found and lower mass limits around 90 — 100 GeV 
were set. Refs. [sj, HH investigated the decay length 
of slepton LSPs assuming trilinear as well as bilinear 
i?-parity violating interactions. Finally, in Ref. [1^, 
the signature of Eqs. ([T|) and ^ was pointed out. But 
in contrast to this work, no signal over background 
analysis was performed. 

The remainder of this paper is organized as follows. 
In Sec.|TT]we review the B3 mSUGRA model and show 
how a in LSP can arise. We present the B3 mSUGRA 
parameter regions with a LSP and propose a set of 
benchmark points for LHC searches. We then classify 
in Sec. IIIII the £r LSP signatures at hadron colliders 



In the B3 mSUGRA model the boundary conditions 
at the GUT scale (A/qut) are described by the six 
parameters [46ll48j 



Afo, Ml/2, Aq, tan^, sgn(^), A. 



(3) 



Here, Mq, Mi/ 2 and Aq are the universal scalar mass, 
the universal gaugino mass and the universal trilinear 
scalar coupling, respectively, tan /? denotes the ratio 
of the two Higgs vacuum expectation values (vevs), 
and sgn{fi) fixes the sign of the bilinear Higgs mixing 
parameter fi. Its magnitude is derived from radiative 
electroweak symmetry breaking [g^I- A is described 
below. 

In B3 mSUGRA, the superpotential is extended by 
the lepton number violating (LNV) terms [68| . 

Wlnv = ^XijkL,LjEk + X-jkUQjDk + n^L^H2, (4) 

which are absent in the MSSM. Here, Li and Qi de- 
note the lepton and quark SU{2) doublet superfields, 
respectively. H2 is the Higgs SU (2) doublet superfield 
which couples to up-type quarks, and Ei and Di de- 
note the lepton and down-type quark SU{2) singlet 
superfields, respectively. i,j,k G {1,2,3} are genera- 
tion indices. Xijk is anti-symmetric in the first two in- 
dices {i O j) and thus denotes nine, A^^j. twenty-seven 
dimensionless couplings. The bilinear lepton number 
violating couplings Hi are three dimensionful param- 
eters, which vanish in B3 mSUGRA at Mqut due 
to a redefinition of the lepton and Higgs superfields 
[46j . However, they are generated at lower scales via 
RGE running with interesting phenomenological con- 
sequences for neutrino masses [2^ [s^l ■ 

In the B3 mSUGRA model, we assume that exactly 
one of the thirty-six dimensionless couplings in Eq. Q 
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is iioii-zcro and positive at the GUT scale^. The pa- 
rameter A in Eq. Q refers to this choice, i.e. 

Ae{A,,fe,A^J, z,j,fc = 1,2,3. (5) 

Given one couphng at the GUT scale, other cou- 
pUngs that violate only the same lepton number 
are generated at the weak scale, Mz, by the RGEs 
[4i,l53n70|,i7ii. 

B. The Selectron and Smuon LSP 

1. Renormalization Group Evolution of the £r Mass 

In order to understand the dependence of the right- 
handed slepton ^, in, mass at Mz on the boundary 
conditions at Mqut, we have to take a closer look 
at the relevant RGEs, which receive additional con- 
tributions from the LNV terms in Eq. The domi- 
nant one-loop contributions to the running mass of the 
right-handed slepton of generation fc = 1, 2 are [i^ 




+ ^>H,k [i^Lh + i^L)n + (mE')fefe] (6) 

with 

(hEk)ij = Xijk X Ao at Mgut, (7) 

and 

S = Tr[mQ2 _ niL^ - 2nif,^ + mg," + mg^] 

+ (8) 

Here, gi (Mi) is the U{1) gauge coupling (gaugino 
mass) and t — InQ with Q the renormalization scale. 
(hEk)ij is the trilinear scalar soft breaking coupling 
corresponding to Xijk. The bold-faced soft mass pa- 
rameters in Eq. © and Eq. ([SJ are 3x3 matrices 
in flavor space: iiIq and m£ for the left-handed dou- 
blet squarks and sleptons, m^, mg, and mg for the 
singlet up-squarks, down-squarks and sleptons, respec- 
tively. rriHi and mjjj are the scalar Higgs softbreaking 
masses. 

The first two terms on the right-hand side in Eq. ^ 
are proportional to the gauge coupling squared, g^, 



^ On the one hand, bounds on products of two different cou- 
plings are in general much stronger than on single couplings 
[69t . On the other hand, one observes also a large hierarchy 
between the Yukawa couplings within the SM. 

^ We consider only the first two generations of sleptons, i.e. 
Ifl S {eij,A-B}i because a stau LSP can also be obtained 
without (large) _R- parity violating interactions [4^. |48|| . 



LiLjEk 


LSP candidate 


2a bound 


Al21, Ai3i 


ea 


0.020 X (Me-j^/100 GeV) 


A231 


cr 


0.033 X (Mg^/100 GeV) 


A132 




0.020 X (A/ji^/100 GeV) 



TABLE L List of LiLjEk couplings (first column) needed 
to generate a cr- or /iij-LSP (second column). The third 
column gives the most recent experimental bounds [95% 
confidence level (C.L.)], taken from Ref. [t^]. The bounds 
apply at AfcuT. The bounds on A212 and A232 from the 
generation of too large neutrino masses are in general too 
strong to allow for a jXR-USP [4^, although exceptions 
might exist [29j . 



and also present in _R-parity conserving models. The 
sum of these two terms is negative at any scale and thus 
leads to an increase of Ma when running from Mqut 
down to Mz. Here, the main contribution comes from 
the term proportional to the gaugino mass squared. 
Ml, because S is identical to zero at Mqut for uni- 
versal scalar masses. Moreover, the coefficient of the 
Ml term is larger than that of the S term. 

The remaining contributions are proportional to 
■^%k ^^"^ ^&)ij': the latter implies also a propor- 
tionality to X^^f. at Mqut, c/. Eq. ([7]). These terms 
are positive and will therefore reduce Mnk , when go- 
ing from Mq\jt down to Mz. They are new to the 
B3 mSUGRA model compared to i?-parity conserv- 
ing mSUGRA. We can see from Eq. (jB]), that if the 
LNV coupling is roughly of the order of the gauge cou- 
pling gi, i.e. Xijk ^ 0(10"^), these terms contribute 
substantially. Then, the in can be lighter than the 
lightest neutralino, Xij and lightest stau, fi, at Mz, 
leading to a Ir LSP [S^. 

The respective LiLjEk couplings A, which can lead 
to a cr or jlR LSP, are given in Table U with their most 
recent experimental 2tT upper bounds at Mqut f72j . 
Because of its RGE running, A at Mz is roughly 1.5 
times larger than at Mqut 0, IZ^ • 

As an example, in Fig. [U we demonstrate the im- 
pact of a non- vanishing coupling A231 at Mgut on the 
running of the cr mass. Note that we can obtain a Cfl, 
LSP {jiR LSP) with a non-zero coupling A121 or A131 
(A132) at A/gut hi a completely analogous way. We 
employ S0FTSUSY3 . 0.13 [zl Izll for the evolution of 
the RGEs. We have chosen a fairly large absolute value 
of = -1000 GeV (see the discussion in Sec. lIIB2|) . 
The other mSUGRA parameters are Mq = 150 GeV, 
Ml/2 = 500 GeV, tan/3 = 5 and > 0. In the corre- 
sponding i?-parity conserving case (A231IGUT = 0), the 
Xi is the LSP and the n is the next-to LSP (NLSP). 

The cr mass decreases for increasing A231, as de- 
scribed by Eq. Furthermore, the masses of the 
(mainly) left-handed second and third generation slep- 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

A231 at A/qut 

FIG. 1: Masses of the en, xi, ti, T2, Vt, fj-L and at 
Mz as a function of A231 at Mgut- The other mSUGRA 
parameters are Mq = 150 GeV, M1/2 = 500 GeV, Aq = 
-1000 GeV, tan/3 = 5 and p > 0. The yeUow region 
corresponds to the experimentaUy aUowed cr LSP region. 
The gray patterned region is excluded by the upper bound 
on A231 , cf. Table [J 



tons, /ii, f2, and sneutrinos, iJ^, Vt-, decrease'^, since 
these fields couple directly via A231. In contrast, the 
mass of the Xi is not changed, since it does not cou- 
ple to the A231 operator at one loop level. Also the 
impact on the mass of the fi, which is mostly right- 
handed, is small. We therefore obtain in Fig. [T] at 
A23i|gut 0.05 a right-handed selectron as the LSP. 

Because of the experimental upper bound on A231 
(see Table H]) the gray pattered region in Fig. [1] with 
A231IGUT > 0.064 is excluded at 95% C.L.. Note, that 
the valid parameter region with a Sr LSP becomes 
larger once we consider scenarios with heavier spar- 
ticles. Moreover, once we go beyond the mSUGRA 
model and consider non-universal masses, a LSP 
can also be obtained with much smaller LNV violating 
couplings. The collider study that we present in this 
publication also applies to these more general £r LSP 
models, provided that we still have a non- vanishing 
and dominant LiLjEk operator. 

In the following, we investigate which other condi- 
tions at Mgut are vital to obtain a £r LSP within 
B3 mSUGRA. Especially the dependence on the tri- 
linear scalar coupling strength ^0 plays a crucial role. 



However, these (negative) ij-parity violating contributions are 
always smaller than those to the right-handed slepton mass 
[4^. Thus, the left-handed sleptons and sneutrinos cannot 
become the LSP within B3 mSUGRA with XijklcuT 7^ 

El. 



2. Aq Dependence 

According to Eq. ([6]) and Eq. d?]), Aq enters the 
running of Ma via the LNV soft-breaking trilinear 

scalar coupling (hEk)ij. As t = liiQ is decreased, 
the (hEk)ij-term gives a negative contribution to M^k ■ 
Its full contribution is proportional to the integral 
of (hEk)-j- over t, from tmin = ln{Mz) to t^ax = 
In(MGUT)- We now show that a negative ^0 with 
a large magnitude enhances the (negative) A^-fc con- 
tribution to the M^k mass. This discussion is similar 

to the case of a sneutrino LSP [ioj . 

In Fig. [2] we show the running of the trilinear cou- 
pling (hgkjij [Fig. 2(a) and the resulting running for 
(h^k)^^- [Fig. 2(b)| . We assume a non- vanishing cou- 
pling AijfcjGUT ~ 0.1 and a universal gaugino mass 
^'Ii/2 = 1000 GeV. Different lines correspond to dif- 



ferent values of Aq, as indicated in Fig. 2(b) 



The dominant contributions to the RGE of (hEk)y 
are given by [46j 



167r 



2 dO^Ek 



dt 



(hEk)y 

18 



ijk 



r.9? + 3gl 



gfMi + &glM2 



(9) 



Ml and M2 are the U(\) and SU{2) gaugino masses, 
respectively. The running in Eq. ^ is governed by two 
terms with opposite sign; one proportional to (hEk)^ 
and one proportional to Ai,fc . In contrast to the sneu- 
trino LSP case {cf. Ref. |49|) the running is indepen- 
dent of the strong coupling 53 and the gluino mass 
M3. 

According to Eq. ([7]), the sign of the term propor- 
tional to (hEk)y in Eq. © depends on the sign of Aq. 
At Mgut, this term is positive (negative) for negative 
(positive) Aq. Hence, for positive Aq, the term propor- 
tional to {h.Ek.) ij increases (hEk)^^ when we run from 
Mgut to Mz. Note, that the gauge couplings gi and 
52 decrease from Mgut to Mz. 

Assuming A^fc to be positive, the second term is 
always positive and thus decreases (hEk)ij when run- 
ning from A/gut to Mz- The Aij^ coupling increases 
by roughly a factor of 1.5 when we run from Mgut to 
Mz. However, at the same time, the gaugino masses 
Ml and M2 as well as the gauge couplings gi and 52 
decrease. Therefore, this term gets relatively less im- 
portant towards lower scales. 

Now, we can understand the running of {h.-£,^)ij in 



Fig. 2(a) Given a positive Aq (red lines), both terms 
in Eq. (|9]) have opposite signs and thus partly com- 
pensate each other, resulting only in a small change of 
(hEk)ij during the running. Moreover, due to the run- 
ning of the gauge couplings and gaugino masses both 
terms in Eq. ^ decrease when we run from A/gut 
to Mz- In contrast, if we start with a negative Aq 
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FIG. 2: Running of (left) and (hEk)?j (right) fr om Mgut to Mz for different values of Aq given in Fig. |2(b)[ At 

A^GUT, we choose A/1/2 = 1000 GeV and Xijk — 0.1. 



(black lines), both terms give negative contributions 
to the running of {h^k)ij. Still, the magnitude of the 
Xijk term in Eq. ^ decreases. Ho-wever, the contri- 
bution from the term proportional to (hEk)y does not 
necessarily decrease -when running from Af gut to Mz ■ 
Thus, for negative Aq, (h^k)ij decreases -with a large 
slope. 



cal uncertainties of SOFTSUSY [78H80| -which is used to 
calculate the SUSY and Higgs mass spectrum. 

The entire displayed region fulfills the 2a constraints 
on the branching ratio of the decay b ^ [8l| , 



3.03 X 10""* < Bib 57) < 4.07 x 10' 



(10) 



Recall from Eq. that Ml is reduced propor- and the upper limit on the flavor changing neutral 



tional to the integral of (h^k)? over t. Thus, accord- 



ing to Fig. 2(b)[ a negative value of Aq leads to a 
smaller Ma compared to a positive Aq -with the same 
magnitude. 



3. Selectron and Smuon LSP Parameter Space 

In this section, -we present t-wo dimensional B3 
mSUGRA parameter regions -which exhibit a In LSP. 
As -we have seen in Sec. IIIB 11 the running of the ^Ii 
mass is analogous for the first and second generation. 
Therefore, -we only study here the case of a cr LSP 
-with a non- vanishing coupling A231 at Mqut- We can 
obtain the fiyi LSP region by replacing coupling A231 
-with A132. 

Wc give in Fig.|3]the LSP region in the Aq — Mi^2 
plane [Fig. 3(a)| and Mq — tan /3 plane [Fig. |3(b)| for a 
coupling A231IGUT = 0.045. We show the mass differ- 
ence, AM, bet-ween the NLSP and LSP. For the sho-wn 
region a lo-wer bound of 135 GeV on the selectron mass 
is employed to fulfill the bound on A231; cf. Table HI 
The pattered regions are excluded by the LEP bound 
on the light Higgs mass [t^, [t^- Ho-wever, wc have 
reduced this bound by 3 GeV to account for numeri- 



currcnt (FCNC) decay 5° 



6(5° 



M^A* ) 



< 3.6 X 10" 



I.e. 



(11) 



at 90% C.L.. 

However, the parameter points in Fig. [3] cannot ex- 
plain the discrepancy between experiment (using pion 
spectral functions from e~^e~ data) and the SM predic- 
tion of the anomalous magnetic moment of the muon, 

; see Ref. and references therein. There exists a 
en LSP region consistent with the measured value of 
at 2a. But this region is alr eady excluded by Teva- 
tron tri-lepton SUSY searches [SJ]. We note however, 
that the SM prediction is consistent with the experi- 
mental observations at the 2a level, if one uses spec- 
tral functions from t data [s^. We have employed 
micrOMEGAs2.2 [sl to calculate the SUSY contri- 



bution to a^, B{b — >■ 57) and S(i3s 



)• 



We observe in Fig. |3] that the e/j LSP lives in an 
extended region of the B3 mSUGRA parameter space. 
Competing LSP candidates are the lightest stau, fi, 
and the lightest neutralino, Xi- 

In the Ao-Mi/2 plane. Fig. 3(a)[ we find a e^ LSP 
for larger values of M1/2, because A'/i/2 increases the 
mass of the (bino-like) faster than the mass of the 
right-handed sleptons [86|, 32} ■ We can also see that 
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(a) Ao-Mj/2 plane with B3 mSUGRA parameter 
Mo = 90 GeV, tan/S = 4, sgn(/i) = -|- and A231IGUT = 0.045. 
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(b) Mo-tan /3 plane with B3 mSUGRA parameter 
Ml/2 = 450 GeV, Aq = -1250 GeV, sgn(/i) = -|- and 
A231IGUT = 0.045. 



FIG. 3: Mass dilference, AM, between the NLSP and LSP. The LSP candidates are explicitly mentioned. The patterned 
regions correspond to models excluded by the LEP Higgs bound. The white dotted line separates e^-LSP scenarios with 
different mass hierarchies: Mg^j < Mf-^ < M^o (left-hand side) and Mg^ < M^o < Mf-^ (right-hand side). 



a Br LSP is favored by a negative Aq with a large 
niagnitude as discussed in Sec. Ill B 21 In this region of 
parameter space the mass difference between the cr 
LSP and the fi NLSP increases with increasing \Aq\. 
In principle, there can also be a en LSP for a large 
positive ^0, c/. Fig. 2(b) However this configuration 
is disfavored due to a too small light Higgs mass [11]. 
Note, that a negative with a large magnitude nat- 
urally leads to a light top squark, ti, since the top 
Yukawa coupling enters the RGB running of the ti 
mass in a similar way as the Xijk Yukawa coupling 
does for the Ir mass [s^ H^l • This behavior plays an 
important role for the mass reconstruction of the fi, 
cf. Sec. El 

In the A/o~tan/3 plane. Fig. |3(b)[ we find a cr LSP 
for tan/3 < 5 and Mq < 100 GeV. The mass of the n 
decreases with increasing tan /3 while the mass of the 
Br is unaffected by tan /3. Increasing tan /3 increases 
the tau Yukawa coupling and thus its (negative) con- 
tribution to the stau mass from RGE running [Sa . ISTj . 
Furthermore, a larger value of tan /3 usually leads to a 
larger mixing between the left- and right-handed stau. 
Thus, tan /? is a handle for the mass difference of the 
fi and cr. In contrast, Mq increases the masses of all 
the scalar particles like the fi and cr, while the mass 
of the Xi is nearly unaffected by both tan/3 and Mq. 
Therefore, at larger values of Mq we obtain a Xi LSP. 

We find basically two possible mass hierarchies for 
the Br LSP parameter space, indicated by the white 
dotted hue in Fig. |3(a)| and Fig. |3(b)[ Close to the Xi 
LSP region, we observe a Xi NLSP and a fi next-to- 
NLSP (NNLSP), i.e. 

(12) 



However, for most of the parameter space, we have 



(13) 



i.e. the n is the NLSP and the Xi is the NNLSP. For 
some regions with a large mass difference between the 
Xi and the cr LSP, the fiR can even be the NNLSP, 
i.e. we have 



Men < Mf, < Mf,n < M^o, 



(14) 



Mg^ < M^o 



< Mf,. 



where the Xi is the next-to NNLSP (NNNLSP). These 
three mass hierarchies lead to a different collider phe- 
nomenology and will be our guideline in the selection 
of benchmark scenarios. 



C. Benchmark Scenarios 

In order to investigate the LHC phenomenology of 
a Br LSP model in more detail, we select for each 
mass hierarchy, Eq. p^ - p^ . one representative br 
LSP benchmark point. The B3 mSUGRA parameters 
and the masses of the lightest four sparticles of these 
benchmark points, denoted BEl, BE2 and BE3, are 
given in Table HIl All benchmark points exhibit a cou- 
pling A231IGUT = 0-045 (cf. Ta ble [J) a nd fulfill the 
experimental constraints of Sec. Ill B 31 and the con- 
straints from Tevatron tri-lepton SUSY searches [U- 
The supersymmetric mass spectra and branching ra- 
tios are given in Appendix El 

The benchmark points BEl and BE2 both feature 
a fi NLSP. In BEl, the fi is nearly mass degenerate 
with the Br and decays exclusively via A231 into an 
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Bs mSUGRA 




benchmark model 




parameter 




BEl 


BE2 




BE3 


Mo [GeV] 







90 




90 


Ml/2 [GeV] 




475 


460 




450 


Ao [GeV] 




-1250 


-1400 




-1250 


tan P 




5 


4 




4 


sgn(/i) 




+ 


+ 




+ 


A23i|gUT 


0.045 


0.045 


0.045 


light sparticles 












(mass/GeV) 












LSP 




(168.7) 


CR (182.3) 




(182.0) 


NLSP 


fi 


(170.0) 


n (189.0) 


X? 


(184.9) 


NNLSP 




(183.6) 


X? (189.5) 


fi 


(187.2) 


NNNLSP 


xS 


(195.7) 


m (199.0) 


P-R 


(195.9) 



A coupling LSP decay 



LHC signature 



TABLE IL Bs mSUGRA parameter and the masses of the 
four lightest SUSY particles of the Sr LSP benchmark 
points BEl, BE2 and BE3. The complete mass spectra 
and the branching ratios are given in Appendix 1X1 



electron and a muon neutrino. In contrast, in BE2 the 
(mainly right-handed) fi is 7 GeV heavier than the 
cr LSP and thus mainly decays via three-body decays 
into the cr due to larger phase-space. Similarly, the 
p,fi NNLSP in BEl decays via three-body decays into 
the CR or the f i . The three-body decays of the heavier 
supersymmetric sleptons to the cr LSP are new and 
are calculated in Appendix |D] 

In BEl, there is a fairly large mass difference be- 
tween the cr LSP and the Xi NNNLSP of about 
27 GeV. The mass difference between the cr and Xi is 
smaller in BE2 (compared to BEl), i.e. about 7 GeV. 
The NNNLSP is the fiR. Finally, the benchmark point 
BE3 features a Xi NLSP that is 3 GeV heavier than 
the CR LSP. The n is the NNLSP and decays into the 
Xi and a r. The fXR is the NNNLSP. 



III. SELECTRON AND SMUON LSP 
SIGNATURES AT THE LHC 

We now classify the main LHC signatures of sclcc- 
tron and smuon, £r, LSP models, under the simplify- 
ing assumption that each decay chain of heavy SUSY 
particles ends in the LSP and that the LSP decay is 
dominated by only one i?-parity violating operator A, 
cf. Table IH If we assume squark pair production as 
the main sparticle production process^, we obtain as 



* For all of our benchmark points the gluinos are heavier than 
the squarks and dominantly decay to a squark and a quark. 
Thus for gluino pair production we simply obtain two jets 
more per event. 



Al21 



sr 



2e 
I 2m 




e Ut 

T Ue 



2e 
[ 2t 



TABLE III: LHC signatures (right column) for selectron 
and smuon LSP scenarios (second column) assuming one 
dominant LiLjEk operator A (left column) and the SUSY 
cascade of Eq. 1)15^ . 



one of the major cascades 

qq/gg^qq 



(15) 



where g is a squark and j denotes a (parton level) jet. 
The two leptons £ are of the same flavor as the LSP. 
The iR_ LSP will promptly decay via the i?-parity vi- 
olating LiLjEk operator into a charged lepton and a 
neutrino. The resulting collider signatures are clas- 
sified in Table IIIII according to the possible £r LSP 
decays. 

Assuming the SUSY cascade in Eq. (fT5|) . the result- 
ing collider signatures involve two (parton level) jets 
from squark decays, two charged leptons from the neu- 
tralino decay with the same flavor as the LSP, as well 
as additional charged leptons and missing transverse 
energy, ^t, from the LSP decays. Because of the Ma- 
jorana nature of the Xi: every charge combination of 
the two £r LSPs is possible. In what follows, it is im- 
portant to note that the transverse momentum, p^, 
spectrum of the leptons from the decay Xi — ^ ^^R will 
depend on the mass difference between the £r LSP 
and the Xi- For smaller mass differences we get on 
average a smaller lepton px . 

In general, more complicated SUSY production and 
decay processes than Eq. can occur. Fig. Ogives 
an example of (left-handed) squark-gluino production 
followed by two lengthy decay chains. Typically, these 
processes lead to additional final state particles [com- 
pared to Eq. ()15p and Table llllj , most notably 



• additional jets from the production of gluinos 
and their subsequent decays into squarks and 
quarks; cf. the upper decay chain of Fig. 21 

• additional leptons from the decays of heavier 
neutralinos and charginos, which may come from 
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pp ■ 




FIG. 4: Example for squark-gluino production with suc- 
cessive cascade decay into two e/?, LSPs. Tiie ii-parity 
violating decays are marked by red dots. (Xn)* denotes 
a virtual neutralino. Note that R-parity violating decays 
can occur earlier in the chain. In this case the LSP is not 
produced. See tables of the benchmark branching ratios in 
Appendix 



the decay of left-handed squarks, like in the 
lower decay chain of Fig. SI and 

• additional leptons from a Xi decay into a non- 
LSP right-handed slepton £^ (or lightest stau 
fi), e.g. Xi — >■ , followed by the three-body 

decay — >■ £'~^£^£]^ via a virtual neutralino, 
iXn)*' ^^'3 the upper decay chain of Fig. U for an 
example. Here, £ii is the LSP. 

These three-body slepton decays are special to cr and 
jlfi LSP scenarios. The corresponding decay rates are 
calculated in Appendix [P] and are taken into account 
in the following collider analysis. 

The coupling A in £ji LSP scenarios is of similar size 
as the gauge couplings and thus enables i?-parity vio- 
lating decays with a significant branching ratio of spar- 
ticles which are not the LSP. Thus, not every SUSY 
decay chain involves the LSP. Of particular impor- 
tance are the 2-body i?-parity violating decays of the 
fi [l^, especially in the case when a fi NLSP is nearly 
mass degenerate with the £r LSP, like for the bench- 
mark point BEl; cf. Table IVlIl Furthermore, sneutri- 
nos (left-handed charged sleptons) may decay into two 



hard charged leptons (one charged lepton and a neu- 
trino) if they couple directly to the dominant i?-parity 
violating operator. This leads to a sharp sneutrino 
mass peak in the respective dilepton invariant mass 
distribution as we will show in Sec. |Vl From the im- 
parity violating left-handed slepton decays we expect 
large amounts of missing energy from the neutrino. 

The lightest top squark, ti, is in most B3 mSUGRA 
scenarios the lightest squark. Thus, ti pair production 
forms a sizable fraction of all SUSY production pro- 
cesses. The decay of each ti yields at least one &-quark 
(either directly from the decay ii — > xt^ and/or from 
the top quark decay after ti — > Xit)- We therefore 
expect an enhanced 6-quark multiplicity for ti pair 
production. We will use the 6-quark multiplicity in 
Sec. |V]to discriminate these events from other SUSY 
processes. 

To conclude this discussion, as one can see from Ta- 
ble mil we expect multi-lepton final states for cr and 
jlu LSP scenarios at the LHC. One the one hand, we 
obtain charged leptons from the Xi decay into the £ji 
LSP. On the other hand, each LSP decay involves a 
charged lepton. Furthermore, as explained above, also 
non LSPs can decay via the dominant _R-parity violat- 
ing operator into leptons. Therefore, a multi-lepton 
analysis will be the best search strategy for our £r 
LSP scenarios. 

Multi charged lepton final states (especially elec- 
trons and muons) are one of the most promising sig- 
natures to be tested with early LHC data. Electrons 
and muons can be easily identified and the SM back- 
ground for high lepton multiplicities is very low [17 1. 



We therefore investigate in the following the discovery 
potential of cr LSP scenarios with an inclusive three 
lepton search analysis. We will treat electrons and 
muons equally and thus expect similar results for fiR 
LSP scenarios. 



IV. DISCOVERY POTENTIAL AT THE LHC 

In this section, we study the discovery potential of 
Br and fiR LSP models with an inclusive search anal- 
ysis for tri- lepton final states at the LHC. Because 
of the striking multi-leptonic signature of these mod- 
els (see Sec. IIII[) . a discovery might be possible with 
early LHC data. We therefore study the prospects at 
the LHC assuming separately a center-of-mass system 
(cms) energy of 7 TeV and 14 TeV. 



A. Major Backgrounds 

In the following Monte Carlo (MC) study, we con- 
sider SM backgrounds that can produce three or more 
charged leptons (electrons or muons) in the final state 
at the particle level, i.e. after (heavy fiavor) hadron 
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and tau Icptoii decays. For the heavy flavor quarks, we 
consider bottom, b, and charm, c, quarks |83- More- 
over, we expect the SUSY signal events to contain ad- 
ditional energy from hard jets arising from decays of 
the heavier (colored) sparticles. We thus consider the 
following SM processes as the major backgrounds in 
our analysis: 

• Top production. Wc consider top pair produc- 
tion (tt), single-top production associated with 
a W boson (Wt) and top pair production in as- 
sociation with a gauge boson {Wti, Zti). Each 
top quark decays into a W boson and a b quark. 
Leptons may then originate from the W and/or 
b decay. 

• Z + jets, i.e. Z boson production in association 
with one or two (parton level) jets. For the as- 
sociated jet(s) wc consider only c- and 6-quarks. 
We force the Z boson to decay leptonically. 

• -I- jets, i.e. W boson production in association 
with two heavy flavor quarks (c or b) at par- 
ton level. We demand that the W decays into a 
charged lepton and a neutrino. 

• Di-boson {WZ, ZZ) and di-boson + jet {WWj, 
WZj, ZZj) production. For the WZ and ZZ 
background, the gauge bosons are forced to de- 
cay leptonically. For WWj, we consider only the 
heavy flavor quarks c and b for the (parton level) 
jet, while for WZj and ZZj every quark flavor 
is taken into account. 

We have also included the processes, where we have a 
virtual gamma instead of a Z boson. 

For the backgrounds with heavy flavor quarks, we 
demand (at parton level) a minimal transverse mo- 
mentum for the c or 6 quarks oi px > 10 GeV cor- 
responding to our object selection cut for the lep- 
tons, cf. Sec. lIVBl Tabic HV] gives an overview of the 
background samples used in our analysis. In princi- 
ple, QCD production of four heavy flavor quarks, like 
bbbb production, can also produce three lepton events. 
However, these backgrounds are negligible compared 
to the other backgrounds in Table IIVI because the 
probability of obtaining three isolated leptons from 
heavy flavor decay is too low [s^ . 

B. Monte Carlo Simulation and Object Selection 

The tt, Zti, Zee and Zbb backgrounds are simu- 
lated with Herwig6.510 [90l-l92t. For the other SM 
processes we employ MadGraph4 .4.30 [9^ for the gen- 
eration of the hard process which is then fed into 
Herwig. The employed MC generators are listed in 
Table [TVl Wc also give the leading-order (LO) cross 



section and the number of simulated events for each 
background sample for both cms energies. The cross 
sections are taken from Herwig (for the tt, Ztt, Zee 
and Zbb backgrounds) or MadGraph (else). We only 
consider the leading-order cross sections for the signal 
and background samples. We note however that the 
next-to-leading-order (NLO) corrections can be large, 
see e.g. Refs. [93 - l98| . and should be included in a 
more dedicated analysis. Furthermore, our simulation 
does not account for detector effects, i.e. we neglect 
backgrounds with leptons faked by jets or photons. 
However, we expect these backgrounds to be small, 
because the fake rate for electrons and muons is quite 
low fl^. 

The SUSY mass spectra were calculated with 
S0FTSUSY3 . . 13 [zl llal • The SOFTSUSY output was 
then fed into ISAWIGl . 200 and ISAJET7 . 64 [g^ in or- 
der to calculate the decay widths of the SUSY parti- 
cles. We added the missing three-body slepton decays 
i'j^ I'dn and fi tUr to the ISAJET code; see 
Appendix [D] for the calculation and a discussion of 
these new slepton decays. The signal processes, i.e. 
pair production of all SUSY particles, were simulated 
with Herwig6 .510. 

We give in Table |V] leading order cross sections for 
sparticle pair production at the LHC for the three 
benchmark scenarios BEl, BE2 and BE3, cf. Tabic HIl 
We separately assume cms energies of ^/s = 7 TeV and 
•y/s — 14 TeV. Wc present the cross sections for the 
signal (last row), i.e. pair production of all sparticles, 
and for three of its subprocesses: The production of 
sparton pairs (second row) , where we consider squarks 
and gluinos as spartons, slepton pair production (third 
row) and the production of electrowcak (EW) gaugino 
pairs or an EW gaugino in association with a squark or 
giuino (fourth row). For all benchmark points, spar- 
ton pair production is the dominant SUSY production 
process. Therefore, the majority of the SUSY events 
will fulflU our signature expectations including at least 
two hard jets. c/. Sec. IHII 

For the rec onstructi on of jets, we employ 
Fast Jet 2.4.1 [loS Il0l| using the kt-algorithm with 
cone radius Ai? = 0.4. Here Ai? = ^TA^F+TA^, 
where rj {(f) is the pseudorapidity (azimuthal angle). 
We only select jets and leptons (i.e. electrons and 
muons) if jryj < 2.5 and if their transverse momen- 
tum is larger than 10 GeV. In addition, leptons arc 
rejected, if the total transverse momentum of all par- 
ticles within a cone of Ai? < 0.2 around the lepton 
three-momentum axis exceeds 1 GeV. 



C. Kinematic Distributions 

In this section we discuss kinematic distributions 
for the benchmark points of Table HIl and motivate our 
cuts of Sec. IIVDI The distributions correspond to our 
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Sample 


Sub- sample 


LO cross 


section [pb] 


Simulated events 


Generator 










7 TeV 


14 TeV 


7 TeV 


14 TeV 








top 


tt 


86.7 


460 


200 000 


5 000 000 


Herwig 








Wt 


10.2 


60.7 


100 000 


1 200 000 


MadGraph 


+ 


Herwig 




wti 


0.14 


0.52 


10 000 


10 000 


MadGraph 


+ 


Herwig 




Ztt 


0.066 


0.43 


10 000 


10 000 


Herwig 






Z + jets 


Zee 


49.5 


187 


100 000 


2 000 000 


Herwig 








Zbb 


44.6 


171 


100 000 


2 000 000 


Herwig 








Z{^ £+£-)+ j {j=e,b) 


59.6 


203 


180 000 


3 700 000 


MadGraph 


+ 


Herwig 


W + jets 


lu)+j3 {j^e,b) 


38.2 


95.2 


135 000 


1 400 000 


MadGraph 


+ 


Herwig 


di-boson 


WZ — >• leptons 


0.20 


0.40 


100 000 


100 000 


MadGraph 


+ 


Herwig 




ZZ — >■ leptons 


0.03 


0.06 


22 000 


75 000 


MadGraph 


+ 


Herwig 




WW + j (j = e,b) 


10.9 


64.0 


120 000 


1 000 000 


MadGraph 


+ 


Herwig 




WZ + j (j = all flavors) 


7.0 


25.0 


77 000 


100 000 


MadGraph 


+ 


Herwig 




ZZ + j {j = all flavors) 


3.2 


10.2 


16 000 


280 000 


MadGraph 


+ 


Herwig 



TABLE IV: SM background MC samples (first and second column) used for our analysis. The third and fourth (fifth 
and sixth) column shows the leading-order cross section (number of simulated events) for pp collisions at a cms energy 
of ^/s = 7 TeV and ^/s = 14 TeV, respectively. For the event simulation we employ the MC generator listed in the last 
column. We also have included the processes, where we have a virtual gamma instead oi a, Z boson. 



Production process 


Cross section [fb] at — 


7 TeV 


Cross section [fb] at ^/s = 


14 TeV 




BEl 


BE2 


BE3 


BEl 


BE2 


BE3 


pp — >■ sparton pairs 


86.7 


152 


139 


1970 


2770 


2760 


pp — >■ slepton pairs 


24.0 


19.9 


21.1 


96.7 


83.9 


88.1 


pp ~> gaugino pairs, 

gaugino-f sparton 


32.2 


38.6 


43.3 


224 


259 


284 


pp — >■ sparticle pairs 


143 


210 


203 


2290 


3110 


3130 



TABLE V: Total LO cross section (in fb) for the benchmark scenarios BEl, BE2 and BE3 for pair production of all SUSY 
particles (last row) and three of its subprocesses: sparton (i.e. squark and gluino) pair production (second row), slepton 
pair production (third row) and electroweak (EW) gaugino pair or EW gaugino plus sparton production (fourth row). 
We separately assume pp collisions at cms energies of = 7 TeV and -^s = 14 TeV. The cross sections are calculated 
with Herwig. We have simulated ^ 15 000 250 000) SUSY events for the 7 TeV (14 TeV) MC signal sample. 



7 TeV event sample and are normalized to one. 

The pt distribution of all electrons [muons] after 
object selection {cf. t he l ast pa ragraph of Sec. IIVBP 
is shown in Fig. [S^a)] [Fig. [5(b)] for the Bg mSUGRA 
benchmark models BEl, BE2 and BE3. In all sce- 
narios, the electrons mostly stem from the neutralino 
decay Xi grc, while many of the muons come from 
the LSP decay iiVt, cf. Appendix 1X1 



We observe in Fig. 5(a) that BEl leads to the in av- 
erage hardest electrons. In this scenario, the mass dif- 
ference between the x? (decaying often via x^ ~^ eci?) 
and the e_R LSP is about 27 GeV and thus quite large 
(compared to the other benchmark points). Further- 
more, the f 1 NLSP decays dominantly via the i?-parity 
violating decay fi — >■ eVf^. A large fraction of the fi 
mass is thus transformed into the 3-momentum of an 
electron. From both sources, we obtain electrons with 
large pt- For example, 81% of all selected electrons 



have p^! > 25 GeV in BEl. 

The situation for BE2 and BE3 is different. Because 
of the smaller mass difference between the Xi and the 
en LSP (compared to BEl), the electrons from Xi de- 
cay are less energetic. For instance, the fraction of 
selected electrons with p^i < 25 GeV is 55% (34%) for 
BE2 (BE3). Furthermore, the electron multiplicity is 
reduced in these scenarios, because many electrons fail 
the lower pT cut (p§! > 10 GeV) of the object selec- 
tion. Due to this, 30% (50%) of all events do not 
contain any selected electron in BE2 (BE3). 



In contrast, the situation for the muons. Fig. 5(b) 
reversed (compared to the electrons) . A large amount 
of the muons arc soft in BEl, whereas BE2 and BE3 
have a harder muon pT spectrum. Note that for BEl, 
a sizable fraction of all muons do not even fulfill the 
object selection requirement of pr > 10 GeV, so that 
34% of all events do not contain any selected muon. 
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0.14- 
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(a)pT distribution of all selected electrons. 



(b)pj' distribution of all selected muons. 



FIG. 5: Electron [Fig. 5(a) and muon [Fig. 5(b) pr distributions at the LHC for the B3 mSUGRA benchmark models 



BEl (red) , BE2 (blue) and BE3 (green) for a cms energy of ^/s = 7 TeV and after object selection cuts. The distributions 
are normahzed to one. 



These muons in BEl stem, for example, from the 3- 
body decays of the p.R into the cr or the fi and are 
in general soft due to decreased phase space, cf. Ta- 
bic IVIII In contrast, the muons in BE2 and BE3 are 
on average much harder, since the majority of these 
muons originate from the cr LSP decay. 

We conclude, that the lepton pT spectrum strongly 
depends on the sparticle mass spectrum. Therefore, 
we desist from making further requirements on the 
lepton PT since this would imply a strong model de- 
pendence in the event selection. We will only require 
at least three charged (and isolated) Icptons as one of 
our cuts in the next section. 

We show in Fig. 6(a) [Fig. |6(b)] the px distribution 
of the [second] hardest jet for the benchmark points 
BEl, BE2 and BE3. For all scenarios, we observe a 
broad peak of the hardest jet pr at around 400 GeV. 
Many of these jets stem from the decays of first and 
second generation squarks into the Xi ; cf. Table IVIII - 
TabUXl 



We find another peak in Fig. 6(a) as well as in 
Fig. |6(b)] at around 100 GeV. These jets stem mainly 
from the t quark decay products from ti — )- txi de- 
cay. The peak is most pronounced in BE2, since here 
we have a light ii mass, M^^ = 448 GeV, and thus 
an enhanced ti pair production cross section. In con- 
trast, the ii mass is about 80 GeV heavier in BEl and 
therefore, the peak is hardly visible in Fig. [51 

For BEl, the pt distribution of the hardest and 
second hardest jet peaks at low values. These soft 
jets stem from initial and final state radiation. They 
appear as the hardest jets in EW gaugino and slep- 
ton pair production which forms a sizable fraction 
(39%) of aU SUSY production processes in BEl, cf. Ta- 
ble |Vl They are less important for BE2 and BE3. 



However, this picture will change for a cms energy of 
y/s = 14 TeV, where sparton pair production is much 
more dominant in BEl. 

Because most events possess at least two jets, we 
demand in the following section at least two jets as 
one of our cuts. Furthermore, we take into account 
that many jets (and some of the leptons) are hard, i.e. 
we demand the visible effective mass to be larger than 
a few 100 GeV; see the next section for details. 



D. Event Selection and Cutflow 

We now develop a set of cuts in order to obtain a sta- 
tistically significant signal and a good signal to (SM) 
background ratio. To motivate the different selection 
steps, we show in Fig. [7] the event distributions that 
correspond to the different cut variables before the re- 
spective cut is applied. We give distributions for the 
three iR LSP benchmark models (BEl, BE2, BE3), for 
the SM background, and, for comparison, for the im- 
parity conserving benchmark model SPSla [5l|^. The 
distributions correspond to an integrated luminosity 
of 1 fb-i at = 7 TcV. 

In Table IVIl we give the number of background and 
signal events after each cut of the analysis. Further- 
more, we provide for each signal benchmark scenario 



^ SPSla has a mass spectrum similar to BEl, BE2, and BE3. 
The main difference lies in the light part of the spectrum, 
where we have in SPSla a stable and invisible Xi LSP. The 
fi is the NLSP. Furthermore, the overall mass scale is a bit 
lower, e.g. the squark and giuino masses are around 500-600 
GeV. 
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FIG. 6: Jet pr distributions at the LHC for the B3 mSUGRA benchmark models BEl (red), BE2 (blue) and BE3 (green) 
for a cms energy of ^/s = 7 TeV and after object selection cuts. The distributions are normalized to one. 



Sample 


Before cuts 


iVlep > 3 


iVjet > 2 


AfoSSF 


> 300 GeV 


top 


97111 ± 197 


14.9 ± 2.2 


13.8 ±2.1 


12.0 ±2.1 


2.1 ±0.8 


Z + jets 


153591 ± 254 


51.9 ±4.3 


16.6 ±2.4 


1.0 ±0.6 


< 1-0 


W + jets 


38219 ± 103 


< 1-0 


< 1-0 


< 1-0 


< 1-0 


di-boson 


21331 ±48 


179.2 ±3.0 


53.5 ±2.0 


2.6 ±0.4 


0.7 ±0.2 


all SM 


310252 ± 341 


264.0 ±5.7 


83.9 ±3.8 


15.6 ±2.2 


2.8 ±0.8 


BEl 


143.1 ± 1.2 


90.5 ±0.9 


79.4 ± 0.9 


68.8 ±0.8 


65.5 ±0.8 


s/Ve 




5.6 


8.7 


17.4 


39.1 


BE2 


210.4 ± 1.5 


92.6 ± 1.0 


81.4 ±0.9 


73.8 ±0.9 


70.4 ± 0.8 


S/^/B 




5.7 


8.9 


18.7 


42.1 


BE3 


202.7 ± 1.4 


61.6 ±0.8 


51.3 ±0.7 


45.2 ±0.7 


43.2 ±0.7 


S/^/B 




3.8 


5.6 


11.4 


25.8 



TABLE VI: Number of SM background and signal events after each step in the event selection corresponding to an 
integrated luminosity of 1 fb^^ at ^/s — 7 TeV. For each signal model (BEl, BE2 and BE3, see Table HI)), we show S/VB 
as significance estimator. The uncertainties correspond to statistical fluctuations. 



the significance estimator S/\/B, where S (B) is the 
number of signal (SM background) events. In general, 
the signal can be defined to be observable if [521 



S > 



5Vb, 5, 0.5B 



(16) 



The requirement S > 0.5B avoids the possibility that 
a small signal on top of a large background could oth- 
erwise be regarded as statistically significant, although 
this would require the background level to be known to 
an excellent precision. In the case of a very low back- 
ground expectation, i? < 1, we still require 5 signal 
events for a discovery. 

As we have seen in Sec. IIIIl we expect an extensive 
number of charged leptons in the final state. However, 
the lepton flavor multiplicity, i.e. the multiplicity of 
electrons and muons, depends strongly on the LSP 
flavor as well as on the dominant A coupling, cf. Ta.- 



ble lllll In addition, as we have seen in the last section, 
the pt spectrum of the leptons is strongly correlated to 
the details of the mass hierarchy. Therefore, in order 
to be as model independent as possible, we simply de- 
mand as our flrst cut three charged leptons (electrons 
or muons) in the flnal state without further require- 
ments on the Pt (beside the object selection cut of 
PT > 10 GeV). 



How useful this cut is, can be seen in Fig. 7(a) 



where we show the lepton multiplicity after object se- 
lection cuts. The distribution for the B3 benchmark 
scenarios peaks around 2-3 leptons, whereas most of 
the SM background events posses less than three elec- 
trons or muons. In principle, by demanding at least 
five charged leptons in the final state, we can already 
get a (nearly) background free event sample. However, 
such a cut would also significantly reduce the number 



13 



no Cut 



s/s=7 TeV 



e 10'' 

« 10''- 
c 
o 

m 10^ 
10 

1 







- BE1 

- BE2 



■ SPS1a 
; SM 



i 9 10 
# Leplons 



(a) Lcpton (e, /j) multiplicity with no cut applied. 



Nlep>2 Cut 



n/s=7 TeV 




10 12 14 16 18 20 
# Jets 

(b) Jet multiplicity after the trilepton cut has been applied. 



NIep+Njet Cut 



\[s=7 TeV 




20 40 60 80 100 120 140 160 180 200 
Invariant Mass of OSSF Leptons [GeV] 



> 

<D 

CD 
o 
o 



NIep+Njet+OSSF Cut 

10F 



\ls=7 TeV 



6- 
4- 
2- 



0, 



• BE1 

• BE2 



-e- SPS1a 
SM 



.-t-=^±-t- 



200 400 600 800 100012001400 1600 1800 2000 

M* [GeV] 



(c) Invariant mass distribution of opposite-sign-samc-flavor (d) Effective visible mass distribution after the trilepton, jet 
(OSSF) lepton pairs after the trilepton and jet multiplicity multiplicity and OSSF mass cut have been applied, 

cut have been applied. 



FIG. 7: Event distributions of several cut variables: Fig. 7(a) lepton multiplicity, Fig. 7(b) jet multiplicity. Fig. |7(c' 



OSSF lepton invariant mass and Fig. 7(d) visible efTective mass for the SM background (gray patterned) and the SUSY 
models BEl (red), BE2 (blue), BE3 (green) and SPSla (white). Note that Fig. 7(a) and 7(b) are given on a logarithmic 
scale. The distributions show the number of events before the event selection cut on the respective variable (see text) 
is applied. They are normalized to an integrated luminosity of 1 fb"'^ at ^/s = 7 TeV. The error bars correspond to 
statistical fluctuations of our MC samples. 



of signal events and is therefore less suitable for an 
analysis of early data. We also observe in Fig. 7(a) 



many more leptons in the i?-parity violating scenarios 
than in SPSla. This is expected, due to the additional 
leptons from the decays of and into the selectron LSP. 

As can be seen in the third column of Table IVTl after 
demanding three leptons, the main SM background 
comes from di-boson events. They account for 68% 
of the background. Furthermore, no + jets events 
survive this cut, indicated by "< 1.0" events in the 
fourth row of Table IVTl At the same time, the number 
of signal events is reduced to 63%, 44% and 30% for 
BEl, BE2 and BE3, respectively. Because of the low 
mass difference between the xl and the bb. LSP in 
BE3, many electrons from Xi decay fail the object 
selection cuts; cf. the discussion of Fig. 5(a) BEl and 



BE2 might already be observable after the first cut, 
I.e. S/v^ > 5. 

Next, we will use the fact that we expect several jets 
from squark and gluino decays; see Sec. IIIII The jet 
multiplicity after demanding three leptons is shown in 
Fig. |7(b)] Because of the weak object selection criteria 
for the jets {pr > 10 GeV) and the small radius for 
the jet reconstruction {AR = 0.4), we observe a high 
jet multiplicity. As discussed in Sec. IIIII we expect at 
least two jets from squark and gluino decays. There- 
fore we demand as our second cut (fourth column of 
Table IVI|) the number of jets to be larger than two, 
i.e. Njct > 2. This cut suppresses roughly two thirds 
of the di-boson backgrounds WZ and ZZ as well as 
of the Z + jets background. However, di-boson pro- 
duction, especially WZ + j, still accounts for most of 
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the background. The number of signal events is only 
reduced by 12%-17%. After this cut, all our bench- 
mark points fulfill the criteria in Eq. and are thus 
observable. 

In order to further reduce the SM backgrounds in- 
volving Z bosons, we construct all possible combina- 
tions of the invariant mass of opposite-sign-same- flavor 
(OSSF) leptons. The distributions (after the three lep- 
ton and A'jet > 2 cut) are shown in Fig. 7(c) As ex- 
pected, the SM background has a large peak at the 
Z boson mass Mz = 91.2 GeV, while the signal dis- 
tribution is mostly flat in that region. Thus as our 
third cut (fifth column of Tabic IVI[) of our event se- 
lection, we reject all events where the invariant mass 
of at least one OSSF lepton pair lies within a 10 GcV 
window around the Z boson mass, i.e. we demand 



MossF ^ [81.2 GeV, 101.2 GeV] . 



(17) 



This cut strongly reduces the Z + jets and di-boson 
backgrounds, leaving ti as the dominant SM back- 
ground. Roughly 90% of the signal events (for all 
benchmark scenarios) survive this cut. The statistical 
significance now lies between 10 and 20 for all bench- 
mark points. 

As we have shown in Sec. IIV C[ our SUSY events 
contain a large amount of energy in the form of high- 
Pt jets and leptons. Thus, we construct the visible^ 
effective mass. 



^"off — 



all 



(18) 



i. e. the scalar sum of the absolute value of the trans- 
verse momenta of the four hardest jets and all selected 
leptons in the event. The visible effective mass distri- 
bution is shown in Fig. |7(d)[ The SM background 
dominates for MJf < 300 GeV, while most of the 
signal events exhibit a visible effective mass above 
300 GeV. This value is slightly higher for the 14 TeV 
dataset. Therefore, we demand as our last cut of our 
event selection (last column of Table IVip 



300 GcV, if^=7TeV, 
400 GeV, ifV^=14TeV. 



(19) 



After this cut, only 2.8 ± 0.8 SM events remain at 
y/s = 7 TeV and an integrated luminosity of 1 fb^^. 
The background is dominated by ti production. The 
signal is nearly unaffected by this cut as can be seen in 
Table IVTl The statistical significance is now roughly as 



® We denote this variable as visible effective mass because it 
does not include the missing transverse energy as in other 
definitions of the effective mass llTtI . 



large as 25 (40) for the benchmark point (s) BE3 (BEl 
and BE2). Furthermore, the signal to background ra- 
tio is now of 0(10). Therefore, systematic uncertain- 
ties of the SM backgrounds are not problematic. A 
signal is clearly visible. 

We observe in Fig. |7(d)| two peaks in the visi- 
ble effective mass distributions for our benchmark 
scenarios. The peak at lower values of M^g con- 
tains mainly events from ti pair production, while 
events from (right-handed) first and second generation 
squark or gluino production build the second peak at 
higher M^^ values. Because of the large mass differ- 
ence between the ti and the other squarks of about 
400 GeV - 500 GeV (depending on the model, see Ta- 
ble |Vll]- Table |IX]), these peaks are clearly separated in 
the visible effective mass. We make use of this fact in 
Sec. IVl when we present a method to reconstruct the 
masses of both the ti and the right-handed first and 
second generation squarks. 

In order to test the flavor sensitivity of our analysis, 
we have applied our cuts to a modified version of the 
benchmark models presented in Table |TT1 Instead of 
A231, we chose A131 (A132) as the dominant i?-parity 
coupling at Mqut to obtain the en (/U_r) as the LSP, 
while leaving the other B3 mSUGRA parameters un- 
changed. The results for the jj-n LSP scenarios are 
in agreement with the original benchmark scenarios 
within statistical fluctuations of the MC samples. 

However, for the LSP scenarios with a domi- 
nant A131 coupling at Mqut, the cut on the invariant 
mass of OSSF leptons rejects more signal events than 
for scenarios with A231. For the modified scenario of 
BEl (BE2), the number of signal events passing the 
-^OSSF cut is reduced by around 15% (3%) compared 
to the original results, cf. Table IVll This difference is 
strongest for BEl-like scenarios, because the endpoint 
of the di-electron invariant mass distribution, where 
one electron comes from the xl decay and the other 
from the br LSP decay, cf. also Eq. (|25ap . coincides 
with the upper value of the Z boson mass window. 
However, this is just a coincidence and a different mass 
spectrum (compared to BEl) with a en LSP and A131 
at Mgut will not have such a suppression. 

We conclude that in most cases, our detailed study 
of e R LSP models with a dominant i?-parity violating 
coupling A231 is representative for all B3 mSUGRA 
models with a en or jlfj LSP. 

To end this subsection, we present in Fig.|8]the miss- 
ing transverse energy, ^t, distribution for the bench- 
mark scenarios, for SPSla and for the combined SM 
backgrounds before any cuts arc applied. In i?-parity 
conserving scenarios like SPSla, the Xi LSP is stable 
and escapes detection leading to large amounts of fx- 
However, for our benchmark points, even though the 
Br LSP decays within the detector, we observe a signif- 
icant amount of missing energy due to the neutrinos 
from the LSP decay. Moreover, the distribution 
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no Cut 
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unaffected. Thus, Mq provides a handle to eontrol the 
mass difference between the xl and the cr (or jlfj) 
LSP. 

the signal cross section (in 
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FIG. 8: Missing transverse energy, fJx, distribution for 
BEl, BE2, BE3, SPSla and the SM background. No cuts 
are applied. The number of events correspond to an inte- 
grated luminosity of 1 fb~^ at y/s = 7 TeV. 



for SPSla falls off more rapidly than in the B3 sce- 
narios. This is because the neutrinos are quite hard, 
resulting from a 2-body decay with a large mass differ- 
ence. Thus, B3 scenarios can lead to even more missing 
transverse energy than i?-parity conserving scenarios. 
We have not employed in our analysis, because 
our simple cuts already sufficiently suppress the SM 
background. Furthermore, it is easier to reconstruct 
electrons and muons than missing energy, especially in 
the early stages of the experiments. 



E. Discovery Potential at the LHC 

In this subsection, we extend our previous analysis. 
We perform a two dimensional parameter scan in the 
M1/2-M0 plane around the benchmark point BEl (Ta- 
ble HT]). For each parameter point, we generate 1000 
signal events, i.e. the pair production of all SUSY par- 
ticles. We then apply the same cuts developed in the 
previous section. We estimate the discovery potential 
of B3 mSUGRA models with a br LSP for the early 
LHC run at -^/s = 7 TeV and also give prospects for 
the design energy of y/s = 14 TeV. 

Due to the RGE running, all sparticle masses at the 
weak scale, especially those of the strongly interact- 
ing sparticles, increase with increasing M1/2 jsil. l87|. 
Thus, by varying we can investigate the discov- 

ery potential as a function of the SUSY mass scale. 
Furthermore, as we have seen in the previous two sec- 
tions, the discovery potential is quite sensitive to the 
mass hierarchy of the lighter sparticles and, in partic- 
ular, to the mass difference between the Xi snd the 
e_R LSP. Increasing Mq increases the masses of the 
scalar particles, while the gaugino masses are nearly 



We show in Fig. 9(a) 
pb) for the LHC with = 7 TeV and in Fig. [9(b)] 
the respective signal efficiency, i.e. the fraction of sig- 
nal events that pass our cuts. The results are given 
only for models with a cr LSP, while models with a 
Xi LSP (fi LSP) arc indicated by the striped (check- 
ered) region. The solid gray region (lower left corner 
of Fig. in]) is excluded by the experimental bound on 
the A231 coupling, cf. Tab|T] 

The signal cross section. Fig. 9(a)[ which is domi- 
nated by the production of colored sparticles, clearly 
decreases with increasing M1/2, i.e. with an increas- 
ing SUSY mass scale. For instance, increasing Afi/2 
from 400 GeV to 500 GeV reduces the cross sec- 
tion from 0.6 pb to 0.1 pb, while the right-handed 
squark (gluino) mass increases from around 820 GeV 
(930 GeV) to 1010 GeV (1150 GeV). In contrast, the 
Mq dependence of the signal cross section is negligible, 
over the small range it is varied. 

For the benchmark scenario BEl, we find in 
Fig. |9(b)| a signal efficiency of 46%. Going beyond 
BEl, we observe that the signal efficiency lies between 
30% and 50% for most of the cr LSP parameter space. 
Therefore, our analysis developed in Sec. IIVDI works 
also quite well for a larger set of cr LSP models. 

However, the signal efficiency decreases dramati- 
cally if the mass difference, AM, between the Xi 
and the cr LSP approaches zero. For models with 
AM < 2.5 GeV, the signal efficiency lies just around 
10% - 20%. As described in detail in Sec. HVCl 
the electrons in this parameter region from the decay 
Xi ~^ Ci{e are usually very soft and thus tend to fail 
the minimum px requirement of the object selection, 
i.e. pT > 10 GeV. For models with AM > 10 GeV, 
the signal efficiency becomes more or less insensitive 
to AAf . Note that, if we choose a stronger minimum 
lepton pt requirement in our analysis, the band of low 
signal efficiency will become wider. 

The signal efficiency depends also slightly on M1/2. 
At low values, M1/2 ^ 400 GeV, i.e. for models with 
a light sparticle mass spectrum, more events are re- 
jected by the cut on the visible effective mass. More- 
over, the SM particles from cascade decays and LSP 
decays have in this case on average smaller momenta 
than in scenarios with a heavier mass spectrum, and 
thus may fail to pass the object selection^. The sig- 
nal efficiency is highest for values of A/1/2 between 
450 GeV and 550 GeV and reaches up to 50%. How- 
ever, when going to very large M1/2, the signal effi- 



However, due to our rather weak requirements for jets and 
leptons, this effect does not play a major role. 
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LHC at V« = 7 TcV 



LHC at = 7 TeV 
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Ml/2 [GeV] 

(a) Signal cross section in pb at the LHC at yTs = 7 TcV. 
The white dashed contour lines give the (first and second 
generation) qji mass in GeV. 



350 400 450 500 550 600 650 700 

[GeV] 

(b) Selection efficiency for the signal events at the LHC at 
y/s = 7 TeV. The white dashed contour lines give the mass 
difference, between the x? and the e/j LSP in GeV. 



FIG. 9: Signal cross section (in pb) [Fig. 9(a) and signal efficiency [Fig. 9(b) at the LHC at ^/s = 7 TeV in the 
— Mo plane. The other parameters are those of BEl (Aq = —1250 GeV, tan/3 = 5, sgn{fi) = +, A23i|gut = 0.045). 
The patterned regions correspond to scenarios with either a fi or xi LSP. The solid gray region in the lower left-hand 
corner is excluded by the bound on A231 , cf. Tab [H 




Ml/2 [GeV] A/1/2 [GeV] 

(a) Minimal required integrated luminosity to obtain an (b) Contours for the signal over background ratio, S/B (red 
observable signal, cf. Eq. III6I I. We show contours for solid lines). 

lOOpb-i (blue dashed line), 500 pb'^ ( green dashed line) 
and Ifb^^ (red dashed line), respectively. 



10: Discovery reach at the LHC at ^/s = 7 TeV in the M1/2-M0 plane. The other B3 mSUGRA parameters are 
— 1250 GeV, tan/3 = 5, sgn(^) = + and A23i|gut = 0.045. We give the minimal required integrated luminosity 

S/B, in Fig. 10(b) The patterned regions correspond 



FIG 

^0 = 

for a discovery in Fig 



10(a] 



and the signal to background ratio 
to scenarios with either a fi or Xi LSP. The solid gray region in the lower left-hand corner is excluded by the bound on 
A231 , cf. Table H] Gray dashed contour lines give the en mass (in GeV) as indicated by the labels. 



ciency again decreases. Here, the production of spar- 
ton pairs is suppressed due to their large masses and 
the jet multiplicity is reduced. Less events will then 



pass the A'jot > 2 and cut. For example, for 

Ml/2 = 500 GeV (M1/2 = 700 GeV), sparton pair 
production contributes (only) 58% (24%) to the total 
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LHC at Vi = 14 TeV 



LHC at = 14 TeV 
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Ml/2 [GeV] 

(a) Signal cross section in pb at the LHC at yfs = 14 TeV. 
The white dashed contour lines give the (first and second 
generation) qji mass in GeV. 



350 500 650 800 950 1100 

Ml/2 [GeV] 

(b) Selection efficiency for the signal events at the LHC at 
y/s = 14 TeV. The wrhite dashed contour lines give the mass 
diflference between the and the LSP in GeV. 



FIG. 11: Same as Fig. [9l but for a cms energy of y/s = 14 TeV. 
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(a) Minimal required integrated luminosity to obtain an 
observable signal, c/. Eq. JTHJ. We show contours for 
lOOpb^^ (blue dashed line), lfb~^ (green dashed line) and 
10 fb^^ (red dashed line), respectively. 
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(b) Contours for the signal over background ratio, S/ B (red 
solid lines). 



FIG. 12: Same as Fig. [TOl but for a cms energy of ^ = 14 TeV. 



sparticle pair production cross section. 

We give in Fig. |10(a)| the discovery potential of cr 
LSP scenarios at the LHC with = 7 TeV. The dis- 
covery reach for the integrated luminosities 100 pb~^, 
500 pb~^ and 1 fb~^ is shown. We use Eq. (fTB)) as 
criterion for a discovery. Furthermore, we present in 
Fig. |10(b)| the signal to background ratio, S/B, as a 
measure for the sensitivity on systematic uncertainties 
of the SM background. As shown in the previous sec- 
tion, the SM background is reduced to 2.8 ±0.8 events 



when we employ the cuts of Table I VII 



Fig. 10(a)| suggests that LSP scenarios up to 



Ml 12 < 620 GeV can be discovered with an inte 
grated luminosity of 1 fb^"'^. This corresponds to 
squark masses of 1.2 TeV and ei? LSP masses of 
around 230 GeV. For these models, we have a sig- 
nal over background ratio of S/B « 3 and thus, sys- 
tematic uncertainties of the SM background are not 
problematic. Furthermore, we see that BEl (A/1/2 = 
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475 GeV, Mq ~ GeV) can already be discovered with 
< 100 pb~^ of data. We also see in Fig. [TU]tliat sce- 
narios with a small mass difference between the Xi smd 
the Br LSP are more difficult to discover as expected 
from Fig. [9(b)| 

We now discuss the prospects of a discovery at the 
LHC at = 14 TeV. I n Fig.lll(a)[ wc give the signal 
cross section and in Fig. 11(b) [ the signal efficiency. We 
employ the cuts developed in Sec. IIVDI The cutflow 
at ^/s = 14 TeV for the benchmark scenarios can be 
found in Appendix IB] 

Because of the higher cms energy, the cross section is 
0(10) times larger than for = 7 TcV, cf. Fig.|9(a) 
For instance, at M1/2 = 400 GeV (500 GeV) the sig- 
nal cross section at y/s = 14 TeV is now 7.2 pb~^ 
(1.7 pb~^). Further more, the signal, i.e. sparticle 
pair production, is now always dominated by sparton 
pair production, cf. also Table [V] 

The signal efficiency at \/s = 14 TeV is slightly 
improved compared to ^/s = 7 TcV. Because of 
the enhanced sparton pair production cross section, 
more signal events pass our cut on the jet multiplic- 
ity, iVjet > 2, cf. also Appendix |B] We now obtain 
a signal efficiency of about 51% (compared to 46% at 
y/s = 7 TeV) for the bench mark p oint BEl. Most of 
the parameter points in Fig. 11(b) exhibit a signal ef- 
ficiency in the range of 40% to 60%. For the scenarios 
with low mass difference between the Xi ^-nd the cr 
LSP, AM < 2.5 GcV, the signal efficiency is reduced 
to around 15% - 25%. As for = 7 TcV, the signal 
efficiency decreases at very large values of A/1/2, be- 
cause of the increasing sparton mass and the reduced 
sparton pair production cross section. Here, this ef- 
fect slowly sets in at values M1/2 ^ 1100 GeV, i.e. 
for scenarios with squark and gluino masses around 
2 TeV. However, even at M1/2 = 1100 GcV, sparton 
pair production still forms half of the total signal cross 
section. 



We show in Fig. 12(a) the discovery potential for 
the LHC at = 14 TeV. We give the discovery 
reach for integrated luminosities of 100 pb~^, 1 fb~^ 
and 10 fb~^, respectively. Our cuts of Scc. lIVDl reduce 
the SM background to 64.7 ± 7.2 events for an inte- 
grated luminosity of 10 fb~^; see Table 1X1 We observe 
that scenarios with M1/2 < 1 TeV (1.15 TeV) can be 
discovered with 1 fb~^ (10 fb~^). This corresponds to 
squark masses of around 1.9 TeV (2.2 TeV) and LSP 
masses of roughly 370 GeV (450 GeV) . The respective 
signal over background ratio is 2 (0.6) as can be seen 
in Fig. 12(b) Therefore, systematic uncertainties of 



the SM background estimate are still not problematic 
as long as the SM events can be estimated to a preci- 
sion of 0(10%). This is a reasonable assumption after 
a few years of LHC running. 

We conclude that due to the striking multi-lepton 
signature, the prospects of an early discovery of B3 
mSUGRA with a iji LSP are better than for i?-parity 



D 



FIG. 13: Decay chain assumed for the mass reconstruction. 



conserving mSUGRA models Note that the vast 
reach in M1/2 is also due to the typically light ti which 
has a large production cross section. For instance, at 
= 525 GeV, the ii mass is around 630 GeV and 
thus can still be produced numerously at the LHC at 
^/i = 7 TeV. 

We want to remark that for scenarios with a low 
mass difference between the x? and the ifj LSP, AM < 
2.5 GeV, the search for like-sign di-lcpton final states 
might be a more promising approach (5^, Il02l4l04 |. 
However, a detailed analysis of these search channels 
is beyond the scope of this paper. 



V. MASS RECONSTRUCTION 

We have shown in the previous section that large 
regions of the B3 mSUGRA parameter space with a 
la LSP can already be tested with early LHC data. 
If a discovery has been made, the next step would be 
to try to determine the sparticle mass spectrum. We 
present now a strategy how the sparticle masses can 
be reconstructed. We use the benchmark point BE2 
as an example. We assume an integrated luminosity 
of 100 fb~^ and a cms energy of y/s = 14 TeV in order 
to have enough events for the mass reconstruction. 

The sparticle decay chains cannot be directly re- 
constructed, because the en LSP decays always into 
an invisible neutrino. Thus, we focus on the measure- 
ment of edges and thresholds of invariant mass dis- 
tributions which are a function of the masses of the 
involved SUSY particles. Our strategy is analogous to 
the one, that is widely used to reconstruct the mass 
spectrum in i?-parity conservi ng SUSY where a stable 
X? LSP escapes detection [TtI. [tosjjToo} . 



A. The Basic Idea 

Wc first discuss the general idea of the method. We 
assume the decay chain 



D^Cc^ Bbc Aabc, 



(20) 
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illustrated in Fig.[T31 where the particles D, C, B, and 
A are massive^ and their masses satisfy 

m£) > mc > TUB > niA- (21) 

The particles c, b and a are observable (massless) SM 
particles. Particle A is assumed to be invisible. 

From the 4-momenta of the decay products a, b and 
c, we can form the invariant mass combinations m{,a^ 
iTT-ca, n^cb and nicba- The maximal (denoted "max") 
and minimal (denoted "min" ) endpoints of these dis- 
tributions, 

max ™max ™max ™max j rn""" (99) 

are functions of the (unknown) particle masses in 
Eq. (|2T|) ^. The respective equations are given in Ap- 
pendix [C] [ul • Note that m^'", m™'^ and mj;;^'" are 
always equal to zero. 

A prominent application of this method is the cas- 
cade decay of a left-handed squark in i?-parity con- 
serving SUSY [13, 

QL ^ qxl ^ q^t^^ ^ Qit^JXi- (23) 

Here, the Xi LSP is stable and escapes the detector 
unseen. Note that in i?-parity conserving SUSY, the 
"near" lepton, and "far" lepton, £f, are of the same 
flavor and thus indistinguishable on an event-by-event 
basis. In our scenarios this is not necessarily the case, 
as shown below. 

For our £ji LSP scenarios, we investigate the decay 
chain of a right-handed squark, i. e. 

qR ^ qx°i -> qt^il ^ qi^t'^y. (24) 

The LSP decays into a charged lepton and a neu- 
trino, where the flavor depends on the dominant A 
coupling, c/. Table IIIIl In contrast to the i?-parity 
conserving scenarios, we can actually distinguish the 
near and far lepton if we have A G {A231, A132}. The 
Ir LSP then decays into a charged lepton of different 
flavor from its own. However, we still have to deal 
with combinatorial backgrounds, because we might 
wrongly combine Icptons (and jets) from different cas- 
cades within the same event. 

In the following, we demonstrate our method for 
the en LSP benchmark model BE2 (A23i|gut 7^ 0), 
cf. Table HIl We focus on the case, where the LSP 
decays into a muon (instead of a r) and a neutrino. 
On the one hand, muons are much easier to reconstruct 
than T leptons. On the other hand, muon events have a 



"near" "far" 




Xi Aasi 



FIG. 14: Relevant decay chain of a right-handed squark, 
Eq. (|24p . for the benchmark scenario BC2. The 7l-parity 
violating decay of the LSP via A231 is marked in red. 

higher probability to pass our cuts, cf. Scct. lIVDl The 
relevant cascade decay, Eq. (|24p . is shown in Fig.[T31 It 
yields one jet (at parton level) and two charged leptons 
of different flavor and opposite charge. From these 
objects, we can form the invariant masses me^, m^q, 
m^g and rrie^,. 

In the mass determination, one can leave the mass 
of the neutrino as a free parameter. If one mea- 
sures this parameter consistent with zero, it would be 
an important piece of information towards confirm- 
ing our model. However, once the i?-parity violating 
decay chain of Fig. [TJ] is experimentally verified (or 
assumed), the knowledge of itia = 0, Eq. ([2T|) . sim- 
plifies the equations in Appendix [C] and reduces the 
number of fit parameters by one. The endpoints of 
the invariant mass distributions are then given by 







(25a) 






(25b) 


(m--)2 


=(*4-^4)(^4-^4)A4' 


(25c) 




-=^4-^l«' 


(25d) 




=A4(Af|o~M|J/(2Af|o). 


(25e) 



In BE2 (and more generally in most £fi LSP models 
within B3 mSUGRA), the (mostly right-handed) ti is 
much lighter than the first and second generation qn. 
Therefore, we have typically two distinct squark mass 
scales. This enables a measurement of the ti and (first 
and second generation) qn mass simultaneously, if we 
are able to separate ti and qn production from each 
other^°. This is possible as we show now. 

B. Event Selection 

For the mass reconstruction, we slightly extend our 
cuts developed in Sec. IIVDI for = 14 TeV. Each 



* Particle A does not necessarily need to be massive. In our 
case it is a massless neutrino. 

^ Another variable which can in princ iple be used fo r our sce- 
narios is the Stransverse mass, 171^2 



From now on, qji stands only for right-handed squarks of the 
first and second generation. 
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event has to contain at least one electron and one 
muon with opposite charge. In order to enhance 
the probability of selecting the right muon, i.e. the 
fi from the LSP decay, we require a minimal 
transverse momentum of the muon of p!^ > 25 GeV. 
We then construct all possible opposite-sign-different- 
flavor (OSDF) dilepton invariant masses, me^, of elec- 
trons and muons (with pt^ > 25 GeV). In order to 
reduce combinatorial backgrounds, we subtract the 
dilepton invariant mass distribution of the same-sign- 
different-flavor (SSDF) leptons. Note that this also 
suppresses (R-parity conserving) SUSY background 
processes, where the charges of the selected leptons 
are uncorrelated, because of an intermediate Majorana 
particle, i.e. a neutralino. For example, SUSY decay 
chains involving the cascade pj^ — > /i^Xi ^ /^^^^^^ 
are thus suppressed. 

For the invariant mass distributions containing a 
jet, we design further selection cuts to discriminate 
between ti and qn events. We expect at least two b 
jets in the ti events from the top quark decays. Thus, 
we introduce a simple 5-tagging algorithm in our sim- 
ulation, assuming a 6-tagging efficiency of 60% pTf . 
We demand two tagged b jets for the ii event candi- 
dates while we require that no 6 jet must be present 
for the qn event candidates. Moreover, we use the vis- 
ible effective mass, M^g, as a handle to discriminate 
between ti and qn events, i.e. we impose the cuts 



400 GeV < MJI < 900 GeV for ii events, 

for events. 



900 GeV < MJf 



(26) 



respectively. 

For the construction of invariant mass distributions 
involving quarks, we consider the hardest and second 
hardest jet, ji and j2 in each event, respectively. Due 
to the lighter ii mass, the jets are expected to be some- 
what softer in ti events than in qn events. Therefore, 
for BE2, we choose the following pT selection criteria 
for the jets: 



50 GeV < p^ < 250 GeV 
25 GeV < pi^ 



} 



for ti events, 



(27) 



250 GeV < 
100 GeV < 



I for qR events. 

P t 



The invariant mass distributions rrif.^q, nieq, and 
m^q are now constructed as follows: 

• fnefiq- We take the invariant masses of the op- 
posite sign electron and muon with ji and j2. 
The smaller (larger) value is taken for the edge 
(threshold) distribution. Note that we repeat 
this procedure for all possible combinations of 
electrons and muons. For the threshold distribu- 
tion, we demand in addition the dilepton invari- 



corresponding to the subset of events in which 
the angle between the two leptons (in the cen- 
ter o f ma ss frame of the LSP) is greater than 
tt/2 |105| . In the edge distribution, we require 



TOe^ < m™^^ and employ SSDF subtraction to 
reduce the combinatorial background. 

• rUeq (m^q): We construct the invariant mass of 
all selected electrons (muons with pt^ > 25 GeV) 
with ji and .72 and take the lower value^^. Fur- 
thermore, we require We^x ^ 



-.max 



For these constructions, the dilepton invariant mass 
edge, m™^^, must have already been fitted. We use 
the true value of the dilepton edge, because it can be 
reconstructed to a very high precision, cf. Sec. IV G II 



C. Results 

We now show our results for BE2 for an integrated 
luminosity of lOO fb"^ at ^/i = 14 TeV. We assume, 
that the SM background can be reduced to a negli- 
gible amount ( cf. Appendix |BJ and present only the 
invariant mass distributions for the SUSY sample, i.e. 
pair production of all SUSY particles. We employ the 
cuts described in the last section. We give a rough 
estimate of how accurately the kinematic endpoints 
may be determined and investigate whether the result 
can be biased due to SUSY background processes or 
systematical effects of the event selection. Our dis- 
cussion should be understood as a proof-of-principle 
of the feasibility of the method. It should be followed 
by a detailed experimental study including a detector 
simulation. 



1. Dilepton Invariant Mass 

We show in Fig. [15] the SSDF subtracted dilep- 
ton invariant mass distribution, rrie^. According 
to Eq. (|25a|) . we expect for the cascade decay in 
Fig. [Tlja dilep ton edge at 51.7 GeV [dashed gray line 
in Fig. 15(a)| . The observed edge quite accurately 
matches the expected value and should be observable 
already with a few fb~^. 

For an invariant mass below the dilepton edge, 
the distribution shape slightly deviates from the (ex- 
pected) triangular shape. This is because the br LSP 
can also decay into a neutrino and a r lepton (see Ta- 
ble |lll|, which then decays into a muon and neutrinos. 



ant mass to lie within m™if^ 



/V2< 



Here, we make use of the fact that we can distinguish the near 
and far lepton. However, we have checked that the model- 
independent con struc tion of the variables m^q{„ear//ar) 
proposed in Ref. [105| | leads to similar results. 



21 



OSDF - SSDF 



OSDF-SSDF 
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me, [GeV] 

(a) Dilepton edge, m^^'^ , Eq. I l25al l . The dashed line 
gives the expected value of 51.7 GeV. 
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(b) Mass peak of the tau sneutrino, i/t, due to the 
i?-parity violating decay i>T — > efi. The true mass is 
Mc^ = 309.8 GeV, cf. Table Ivml 



FIG. 15: Dilepton invariant mass distribution, m^^, for the benchmark point BE2. The distributions are same-sign- 
different-fiavor (SSDF) subtracted. The error bars correspond to statistical uncertainties at 100 fb^^. 



In this case, the muon only carries a fraction of the r 
lepton pt and we obtain an on average lower value 
compared to the LSP decay e^, — > fii^r- 

We observe another small edge at about 70 GeV. 
These events stem from the decay of a left-handed 



smuon, I.e. jr^ — > M^Xi 



li'^e^eR, cf. Table Ivml 



The true endpoint is 70.7 GeV. 

Furthermore, as shown in Fig. |15(b) we have a 
sharp peak at 309.8 GeV in the di-lepton invariant 
mass distribution. Here, the mass of the tau sneutrino, 
i^T, is fully reconstructed. It decays via the i?-parity 
violating decay f'r — > e^/i^ with a branching ratio of 
12%; see Table IVIIII Analogously, we also expect a 
mass peak in the er invariant mass distribution from 
the respective muon sneutrino decay. However, the 
observation of this peak requires the reconstruction of 
the T lepton momentum which is beyond the scope of 
this paper. The sneutrino mass peaks are expected to 
be observable with only a few fb~^ of data and are 
thus a smoking gun for our scenarios. 



2. Dilepton plus Jet Invariant Mass 



We show in Fig. (TH] the dilepton plus jet invariant 



^euQ 



, to obtain the kinematic edge 



mass distribution, 

for the qn event [Fig. 16(a) and ii event [Fig. |16(b)| 
selection, cf. Sec. lVBl Recall that we employ different 
selection criteria to obtain the edge and the threshold 
of the TTie^q distribution; see the end of Sec. IV Bl for 
details. 

According to Eq. ()25d|) and Table I Villi we expect 
the edge in Fig. 
[410 GeV]. For t: 



16(a)] [Fig. [16(b)] to lie at 925 GeV 



as can be seen^^ by the dashed gray line in Fig. 16(a) 



In contrast, in the ti event selection the identifica- 
tion of the endpoint [dashed gray line in Fig. |16(b)| 
is more difficult. The observable edge is smeared to 
higher values. On the one hand, cascade decays of 
heavier squarks and gluinos can leak into the ti event 
selection. On the other hand, the distribution flattens 
out as it approaches the nominal endpoint, because 
the jet (from t decay) carries only a fraction of the t 
quark pT. Moreover, the cut imposed on the jet trans- 
verse momentum, px < 250 GeV, Eq. ((27)) . tends to 



reject events at high mefj,q values. Therefore, the end- 
point tends to be smeared. However, the intesection 
of the X-axis with a linear fit on the right fiank of 
Fig. |16(b)] would still provide a quite good estimate of 
the true edge. Such a procedure is also employed for 
the mass reco nstruction of i?-parity conserving models 

In Fig. [ITl we present the mpun threshold- 
[Fig 



distribution for 



we 
the 



17 a 



[Fig. |17(b)] selection. In Fig. [17 



■eiiq 

and 



ti event 
we observe an 

edge slightly below the expected threshold of 181 GeV 
(gray dashed line). This shift towards lower values is 
mainly d ue to final state radiation of the quark from 
qn decay |ll5| , i. e. the reconstructed jet is less ener- 
getic than the original quark. This is not surprising, 
because we use a relatively small radius, Ai? = 0.4, 
for the jet algorithm, cf. Sec. IIVBI 

In general, the m^^q threshold value is set by the 
lightest squark. Therefore, events in Fig. 17(a) with 



values far below the endpoint at 181 GeV usually con- 



re qpi event selection, this is the case 



The endpoint values are usually determin ed by employing 
straight line fits, see e.g. Ref. [rn.ll05| - [l07ll . 
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(a) rriefiq edge-distribution for the qn event selection. 
The dashed hne gives the expected value of 925 GeV, 
cf. Eq. l l^SdT l. 
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(b) mefiq edge-distribution for the ti event selection. 
The dashed line gives the expected value of 410 GeV, 
cf. Eq. ll25dll for the tj . 



FIG. 16: Dilepton plus jet invariant mass distributions, mef_iq, for the kinematic edge for the qn event [Fig. 16(a) and ti 
event [Fig. |16(b)| selection. The distributions are SSDF subtracted. The errors correspond to statistical uncertainties at 
100 fb-\ 
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(a) rrieiiq threshold-distribution for the qn event 
selection. The dashed lino gives the expected value of 
181 GeV, cf. Eq. 



(b) rrieiiq threshold-distribution for the ii event 
selection. The dashed line gives the expected value of 
86 GeV, cf Eq. |[25ct for the h . 



FIG. 17: Dilepton plus jet invariant mass distributions, me 



and ti event [Fig. |17(b) selection. The errors correspond to statistical uncertainties at 100 fb 



for the kinematic threshold for the qr event [Fig. 17(a) 

1 ' 



tain third generation squarks. These events can leak 
into the qr event selection when the b quarks are not 
tagged. 



For the 
served m. 



ti event selection, Fig. 17(b) the ob- 
efj,q threshold matches quite accurately the 
expected value of 86 GeV (gray dashed line). We 
note however, that detector effects, especially jet miss- 
measurements, are expected to smear the thresholds 
and edges. But, this lies beyond the scope of this pa- 
per. 



3. Lepton plus Jet Invariant Masses 



We now discuss the invariant mass distributions 



formed by one charged lepton and a jet. 



''eq 



and 



For these invariant masses, we generally have 
larger SUSY backgrounds (compared to the dilepton 
and dilepton plus jet invariant mass distributions), be- 
cause we cannot employ SSDF subtraction. 

The electron-jet invariant mass distributions, nieq, 

In the qn event selection 
, we observe an edge near the expected end- 
point of 251 GeV (gray dashed line). In contrast, in 



are presented in Fig. 
[Fig. 



1^ 



the ti event selection [Fig. 18(b)| , the endpoint, which 
is expected to lie at 111 GeV, cannot be easily identi- 
fied. 
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(a) rrieq distribution for the qjf event selection. The 
dashed Une gives the expected value of 251 GeV, 
cf. Eq. | [25ct . 
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(b) rrieq distribution for the ti event selection. The 
dashed line gives the expected value, 111 GeV, 
cf. Eq. |[25cll for the ti . 



FIG. 18: Electron plus jet invariant mass distribution, rrieq, for the qa event [Fig. |18(a) and ti event [Fig. |18(b) selection. 
The errors correspond to statistical uncertainties at 100 fb"^. 



The jet used for Fig. 18(b) usually carries only a 
fraction of the t quark momentum reducing the in- 
variant mass. In addition, the ti cascade decay 



L ^HI^ b^+r, h„+ 



t 



bxt 



(28) 



possesses an endpoint at 267 GeV in ra^q which pro- 
duces events beyond the expected endpoint. As a re- 
sult, a measurement of the 111 GeV endpoint will be 
difficult. 

In Fig.[T2]we show the muon-jet invariant mass dis- 
tributions for the qn event [Fig. 19(a)| and ii event 
[Fig. |19(b)| selection. Assuming the cascade de- 



cay of Fig. [TH the m^q distribution. Fig. 19(a) has 



an expected endpoint at 921 GeV, Eq. ()25bl) . We can 
clearly observe an endpoint in Fig. 19(a) However, 
in general it might be slightly underestimated, due to 
final state radiation of the quark from squark decay. 

In the ti event selection, the endpoint is again 
more difficult to observe, cf. Fig. 19(b) For ra^q > 
300 GeV, the distribution approaches the endpoint 
with a very flat slope. Thus, the determination of the 
endpoint requires high statistics. Moreover, we have 
background events beyond the endpoint from heavier 
squark cascade decays or combinations with a jet from 
a decaying gluino. 

We conclude that the standard method that is used 
to reconstruct sparticle masses in i?-parity conserving 
SUSY works also very well for our £ji LSP models, 
where the LSP decays semi- invisibly. We therefore 
expect that most of the SUSY masses in our model can 
be reconstructed with a sirnilar p recision as in i?-parity 
conserving models [T7L [T05l - [l09| . i.e. we expect for our 
benchmark model a relative error of about 10% or less. 
We have not calculated the sparticle masses from the 
kinematic edges, because for a reliable estimate of the 



errors, one has to include detector effects. However, 
this lies beyond the scope of this work. 



VI. SUMMARY AND CONCLUSION 

If i?-parity is violated, new lepton number violating 
interactions can significantly alter the renormalization 
group running of SUSY particle masses if the coupling 
strength is of the order of the gauge couplings. Within 
the framework of the B3 mSUGRA model, we showed 
that a selectron and smuon LSP can arise in large re- 
gions of the SUSY parameter space {cf. Fig. [3]) if a 
non-vanishing lepton number violating coupling Xijk 
with k ~ 1,2 is present at the GUT scale; see Table U 
for a list of all allowed couplings. 

The selectron or smuon LSP decays mainly into a 
charged lepton and a neutrino. Additional charged 
leptons are usually produced via cascade decays of 
heavier sparticles into the LSP. Keeping in mind that 
sparticles at the LHC are mostly produced in pairs, 
we end up with roughly four charged leptons in each 
event at parton-level. Furthermore, two or more jets 
are expected from decays of strongly interacting SUSY 
particles. Table Hill gives an overview of the expected 
LHC signatures. 

Based on this, we have developed in Sec. llVI a ded- 
icated trilepton search for our SUSY scenarios. We 
found that demanding three charged leptons and two 
jets in the final state as well as employing a Z-veto 
and a lower cut on the visible effective mass is suf- 
ficient to obtain a good signal to background ratio. 
For example, for an integrated luminosity of 1 fb^-'^ 
at = 7 TeV, only approximately three SM events 
survive whereas the number of SUSY events passing 
our cuts can be of 0(10 - 100), cf Table Ell 



24 



220 
200 
180 
160 
140 
120 
100 
80 
60 
40 
20 
0, 



q events 



t, events 



200 



400 



600 



I , ,'«i..-,.......L.j.u 

800 1000 
m„„ [GeV] 



(a)mp5 distribution for the gjj event selection. The 
dashed hne gives the expected value of 921 GeV, 
cf. Eq. l l25bt . 
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(b)mfiq distribution for the ti event selection. The 
dashed line gives the expected value of 406 GeV, 
cf. Eq. ll25bll for the tj . 



FIG. 19: Muon plus jet invariant mass distribution, m^q, for the qr event [Fig. 1 19(a) and ii event [Fig. |19(b) selection. 
The errors correspond to statistical uncertainties at 100 fb~^. 



We found within the B3 mSUGRA model that sce- 
narios with squark (selection or smuon LSP) masses 
up to 1.2 TcV (230 GeV) can be discovered with an 
integrated luminosity of 1 fb~^ at y/s = 7 TeV, thus 
exceeding the discovery reach of i?-parity conserving 
models. Our scenarios are therefore well suited for an 
analysis with early LHC data. Going to a cms energy 
of =14 TeV and assuming an integrated luminos- 
ity of 10 ih~\ allows a discovery of 2.2 TeV (450 GeV) 
squarks (sclcctron and smuon LSPs). 

After a discovery has been made, a next step would 
be the reconstruction of the SUSY mass spectrum. 
Unfortunately, although the LSPs decay, a direct mass 
reconstruction is often not possible (sec Fig. |15(b)| 
for an exception), because (invisible) neutrinos are al- 
ways part of the LSP decays. We therefore proposed 
in Sec. |V] a method relying on the measurement of 
kinematic edges of invariant mass distributions. This 
method is analogous to the one usually used for im- 
parity conserving models, although different SUSY 
particles arc involved in the decay chain. For example, 
the neutrino from the LSP decay in our models plays 
the role of the lightest neutralino in i?-parity conserv- 
ing models. We also showed that decay chains from 
heavier (first and second generation) squarks can be 
distinguished from those of the lighter (third genera- 
tion) top-squarks. Therefore, a measurement of both 
squark mass scales is possible. 
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Appendix A: Properties of the Benchmark Models 



We show in Tables IVIIi IVIIII and IIXI the mass spec- 
tra and the dominant decays of the supersymmetric 
particles of the benchmark points BEl, BE2 and BE3, 
respectively; see Table for a definition. Sparticle 
masses, that are reduced by more than 5 GeV (com- 
pared to the i?-parity conserving case) and i?-parity 
violating decays are marked in bold-face. Note that 
only the masses of those sparticlcs, which couple di- 
rectly to the LiLjEk operator, are significantly re- 
duced, cf. Sec. IIIB II Therefore, in our benchmark 



models (A231IGUT 0) only the e^,, /2l, i>^, T2 and Vr 
are affected. 

These sparticles then also exhibit i?-parity violat- 
ing decays to SM particles via A231. For the Vt- this 
can lead to a striking peak in the electron-muon in- 
variant mass distribution; cf. Fig. 15(b) In addition, 
the fi can also decay via the A231 coupling, because 
of its (small) left-handed component. This happens 
in particular in scenarios, where the fi is the NLSP 
and its mass is close to the LSP mass (as in BEl, Ta- 
ble |VlIl), j-e. the i?-parity conserving decay into the 
LSP is phase-space suppressed. The LSP can only 
decay via i?-parity violating interactions: cr (ivt 
and Br Tf^. 

Common to all benchmark points is a rather light ti 
(compared to the other squarks). For all benchmark 
points, the ii mass is around 450 GeV-550 GeV and 
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Mass [GoV] 


Channel BR 


Channel BR 


168.7 


/^"i/t 50% 


T^iy,, 50% 


fj" 170.0 


e^p^ 100% 




jX~^ 183.6 


e^e~/j^ 34.6% 
f + T-fi- 20.4% 


e^e+ 28.3% 
f-r + fi- 16.7% 


X? 195.7 


e~e+ 23.8% 
fj"T+ 21.0% 
/Ip/i+ 5.1% 


e+e" 23.8% 
f^T" 21.0% 
ApM^ 5-1% 


£>^ 306.5 


Xi'i^T 60.2% 
e";i+ 11.4% 


W+f- 28.4% 


£>j, 309.4 


X>f, 84.4% 


e-T+ 15.6% 


313.5 


xVe 100% 




Tj" 318.4 


x'It- 59.0% 
Z"?- 14.1% 


H"f- 16.5% 
e"P„ 10.4% 


/I7 318.7 


x'lf^- 84.1% 


e^Px 15.9% 




Xi*- iuu/0 




xl 372.0 


5^1Jr 10.0% 

P^i/j, 9.2% 

5el^e 8.1% 

f~T+ 7.1% 
e-e+ 6.2% 
f-r+ 1.6% 


PxPx 10.0% 
i^^P^ 9.2% 
PePe 8.1% 
r/+ 7 90^ 
f + T-- 7.1% 
e+e" 6.2% 
f + T-- 1.6% 


xr 372.0 


S^r" 20.6% 
v^e- 16.8% 
Tj^Pt 13.7% 
r~ u 3 1 % 


19.0% 
A^Pf, 13.9% 
e£Pe 12.0% 


ti 531.1 


X?t 62.2% 


X+& 37.8% 


61 847.3 


W^ti 71.5% 
v^/) 10 4% 

X2" J-U.^/O 


Xi^t 17.5% 


X3 898.0 


tit 19.7% 
XiW+ 18.4% 
v^Z" 16 5% 
X^//" 1.2% 


tit 19.7% 

X+W" 18.4% 

v^z" 4 8% 
X 1 ^ ^ . /o 


X2 906.0 


i'lh 47.6% 

-v~ 1 ^ 4% 
X 1 ±0.^/0 

xJW^" 4.2% 


X2l^" 15.9% 

^~ 1 4 fi% 

Xl ^ -L4:.U/0 


-v*^ 906 4 
X4 JUD.^ 


f 2Q fi% 
v^W"*" 12 1% 

x^i/° 10.3% 


-y+W^ 12 1% 

xjff" 2.9% 


t2 919.4 


z^fi 49.1% 

X+b 17.3% 
xjf 1.5% 


H"ii 24.6% 
fSf 7 6% 


62 959.5 


X?& 67.0% 
Xr* 2.1% 


W^fi 28.9% 
X°& 1.2% 


rfn (sr) 962.3 


x;<i(s) 100% 




mr(cb) 965 


x;«(c) 100% 




Mi,(£t) 1001.8 


X+d(s) 65.9% 
X;«(c) 1.2% 


X2«(c) 32.9% 


cii,(st) 1004.7 


X~«(c) 65.5% 
xVls) 1-7% 


X°d(s) 32.8% 




f-, f 9(1 Q% 

&1& 8.5% 

&2& 2.9% 
dRd{sRs) 2.7% 
Ufiu{ciic) 2.6% 
t2t 1.6% 
mi,m(ci,5) 1.4% 
dLd(sLs) 1.3% 


1; -[^ [■ iiu . y /o 

bib 8.5% 
636 2.9% 
(i^d(s^s) 2.7% 
iJi«(c^c) 2.6% 
tjt 1.6% 
i2"(c2c) 1.4% 
J^dCsis) 1.3% 



TABLE VII: Branching ratios (BRs) and sparticle masses 
for the benchmark scenario BEl; see TableUll BRs smaller 
than 1% are neglected. iZ-parity violating decays and 
masses which are reduced by more than 5 GeV (compared 
to the i?-parity conserving case) are shown in bold-face. 



Mass [GeV] 


Channel 


BR 


Channel 


BR 


R 


182.3 




50% 




50% 


T-,~ 

1 


189.0 


e+e~ r~ 


50.2% 


e-e+T- 


49.5% 


X? 

A. 1 


189.5 


e T-, 
'-7?'- 


50% 


e+e- 


50% 


U >-. 
H-R 


199.0 


X?M" 


100% 






309.8 




71.0% 
12.0% 


w+f- 


17.0% 


fi 


312.0 


xVm 


85.8% 


e-T+ 


14.2% 




317.0 


X?l'e 


100% 






320 8 


Xjr- 


69.9% 


e i/j. 


11.3% 








10.2% 


Z 


8.6% 




320 8 


AIM 


85.2% 




14.8% 




325.7 


Xie 


100% 




X2 


360 1 




1 u . /o 




10.5% 








9.7% 




9.7% 








7.9% 


Of, Pe 


7.9% 






^2 ^ 


7.0% 

D.O /O 


r^r- 


7.0% 
6.8% 






e^ e+ 


5.4% 


e+e 


5.4% 








2.0% 


f + r- 


2.0% 






Xl 


1.3% 






~_ 

Xl 


360.2 




21.7% 




19.9% 








16.3% 




13.4% 








13.2% 


e-p. 


10.5% 






'1 " T 


3.8% 


X?M/- 


1.3% 




448 3 


x?t 

Xl 


71.9% 


X^ 


28.1% 


L"! 


809 1 


VK~ fi 
\% 

A 2 


78.1% 
8.1% 


Xl t 


13.3% 




887 


z"ti 


52.7% 


H°ii 


25.9% 






v+6 


14.1% 


x%t 


6.1% 






A,l " 


1.2% 






X? 


936.7 


tit 


26.0% 


tit 


26.0% 






Xi;VF+ 


14.6% 


xtw- 


14.6% 






A2 


13.3% 


X°iZ° 


3.8% 




937.7 


Al " 


67.9% 


W'ti 


26.0% 






Al 


2.7% 




1.5% 


dn. (sr) 


939.8 


Xl "l,'^,' 


100% 




^R (Cr) 


942 9 


x!u(c) 


100% 




X2 


944.5 




55.6% 


xlw- 


13.5% 








13.1% 


XiH° 


12.4% 






Xl 


3.4% 






X4 


945. 1 


t\i 


33.3% 


ilt 


33.3% 






Al 


10.0% 


xtw~ 


10.0% 






X2-H"° 


8.5% 


xlH° 


2.4% 




977.6 


\td(s) 

Al ""V-/ 


65.9% 
1.2% 


Xlu{c) 


32.9% 


O.L ) 


980 4 


Xl '^^(c) 

Xl"'!,'^,' 


65.6% 
1.6% 


Xld[s) 


32.8% 


9 


1063.1 


i\ t 


23 0% 


ilt 


23.0% 






bib 


8.7% 


bib 


8.7% 






626 


2.4% 


Kb 


2.4% 






dRd(sRs) 


2.4% 


d'j^d{s*j^s) 


2.4% 






uru{crc) 


2.3% 


Uru{c1^c) 


2.3% 






t2f 


2.0% 


ilt 


2.0% 






"l,fi(ci,c) 


1.2% 




1.2% 






dLd(sLs) 


1.1% 




1.1% 



TABLE VIII: Same as Table |VlIl but for the benchmark 
scenario BE2. 
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Mass [GeV] 


Channel BR 


Channel BR 


g;; 182.0 


fJ-^U-r 50% 


r^iy^ 50% 


x" 184.9 


ej,e+ 50% 


e+e" 50% 


fi 187.2 


X?!-- 64.5% 


e"P^ 35.5% 


u~ 19.'5.9 


xiV" 100% 




304.3 


Xji^x 73.6% 
e~/i+ 12.2% 


W+f~ 14.2% 


£>j, 306.2 


X>f> 86.0% 


e~T+ 14.0% 


£>e 310.4 


Xi'i^e 100% 




/t^ 315.2 


XiV" 85-2% 


e"pT 14.8% 


f~ 315.3 


X^r- 72.5% 
H°f- 8.5% 


e-p^ 11.7% 
ZOf- 7.3% 


319.3 


xje- 100% 




X2 351.2 


Z-'x i'T- 1 U . /O 

S^Pj, 9.7% 

5e!^e 8.1% 
6.8% 

T2 T D.D/0 

ffT+ 2.0% 
xjff" 1.6% 


I^x J^T IV. /o 
8.1% 

jl^fi- 6.8% 

T D . D /O 
e e . 'rt /o 

f + r- 2.0% 


351.2 


S^T" 21.7% 
See" 16.8% 
Tj^Pt 12.7% 
fj^P^ 3.8% 


20.1% 
A^Pf, 13.0% 
e^Pe 10.4% 
xJW- 1.6% 


ti 481.7 


X?* 62.1% 


X+b 37.9% 


61 805.4 


IV^ti 73.9% 
X 2 y . / /o 


Xi^t 15.9% 


t2 881.7 


Z"*! 51.3% 
X+b 16.1% 
-v^f 1 4% 


if"ti 24.2% 
xlt 7.0% 


X3 884.0 


tit 22.1% 

XrW+ 17.0% 

v° Z° 15 4% 
X2-^ J-O.^/o 

X^ff" 1.0% 


fit 22.1% 
xtW- 17.0% 
v^Z" 4 5% 

Al^ '±.%J/0 


v" 892 4 


Xl ^ i'i.U/O 

xjw 3.8% 


W 1 5 1 % 
ff" 13 7% 


xl 893.1 


tit 31.4% 
X^W"*" 11.1% 
X°ff" 9.4% 


fit 31.4% 
XiW~ 11.1% 
x'lH° 2.7% 


62 919.3 


Xib 70.5% 
xr* 1-8% 


W~ii 26.0% 
X2b 1-1% 


dR{sR) 921.1 


x;d(s) 100% 




MH 923.8 


Xiu{c) 100% 




UL 957.9 


X+rf(s) 65.9% 
X'ln{c) 1.1% 


X2m(c) 33.0% 


rft(sL) 961.0 


xr"(c) 65.5% 

Xl ±. ( /o 


X^(i(s) 32.7% 


g 1041.8 


fit 22.6% 
bi& 8.9% 
Sab 2.7% 
dRdisRs) 2.6% 
uru{crc) 2.5% 
mlm(clc) 1.3% 
dLd{sLs) 1.2% 


fit 22.6% 
b*b 8.9% 
bjb 2.7% 
d*j^d{s*Rs) 2.6% 
ii^u(c^c) 2.5% 
^^^(e^c) 1.3% 
1.2% 



TABLE IX: Same as Table |VlIl but for the benchmark 
scenario BE3. 



the other squark masses are in the range of 800 GeV- 
1 TeV. Because of the large top Yukawa coupling, the 
stop mass receives large negative contributions from 
RGE running, especially for a negative Aq with a large 
magnitude [8a, [S^; see Sec. Ill B 21 for a similar case. 
Furthermore, the light stop mass is reduced by large 
mixing between the left- and right-handed states. As 
one can see in Tables IVIIl I Villi and IIXI the (mainly 
right-handed) ti dominantly decays into the (bino- 
like) Xl s-nd a top quark, while the decay into the 
(wino-like) lightest chargino, xt ^ subdominant. 

The in, flji, fi and Xi always form the lightest four 
sparticles in B3 mSUGRA models with a en or fin 
LSP. The scenario BEl, Tablc|VlTl exhibits a n NLSP 
that is nearly degenerate in mass with the en LSP. 
Thus, it undergoes the i?-parity violating decay fi — >■ 
eVfi, yielding high-pr electrons. The jln is the NNLSP 
and decays into the en or the fi via 3-body decays 
producing in general two low-pr charged leptons due 
to the reduced phase space. We calculate and discuss 
these decays in detail in Appendix |DJ The Xi is the 
NNNLSP. Besides the decay into the en LSP and an 
electron (47.6%), it also decays to a sizable fraction 
(42.0%) into the fi NLSP and a r lepton. 

The benchmark scenario BE2, Table IVTIIl also has 
a fi NLSP. However, the fi NLSP is nearly mass de- 
generate with the Xl NNLSP. Therefore, it decays ex- 
clusively via 3-body decays into the en LSP, yielding 
a low-pT tau lepton and an electron; cf. Appendix [Dl 
The Xl NNLSP always decays into the en LSP and an 
electron. 

In contrast to BEl and BE2, the NLSP in BE3, 
Table IIXI is the Xi which is roughly 3 GcV heavier 
than the en LSP. Therefore, the electrons from the 
Xl decay into the LSP arc very soft. We have a fi 
NNLSP, which decays i?-parity conserving into the x? 
and a tau as well as via i?-parity violating decays into 
an electron and a neutrino. In both BE2 and BE3, the 
jln is the NNNLSP and decays exclusively into the Xi 
and a muon. 

The remaining sparticle mass spectra and decays 
look very similar to those of i?-parity conserving 
mSUGRA [53. 



Appendix B: Cut-Flow for = 14 TeV 

We present in Table 1X1 the cut flow of the signal and 
SM background events for an integrated luminosity of 
10 fb"^ at = 14 TeV. Ahhough the benchmark 
scenarios BEl, BE2 and BE3 (see Tabic HI]) are already 
observable with very early LHC data, cf. Scc. lIVDl we 
provide their expected event yields reference in 
order to compare the signal efficiencies at -^/i = 7 TeV 
and = 14 TeV. 

We apply the inclusive three-lepton analysis devel- 
oped in Sec. IIVDI After the three lepton requirement 



27 







N, > 

-fVlcp ^ 








top 


/coi c _|_ 9\ 1 n3 


000 zn z 1 




■^07 -1- 1 Q 


Q -1- 7 n 


Z + jets 


(5601 ± 2) • 10 


1980 ± 41 


571 ± 22 


48.7 ± 6.4 


2.6 ± 1.5 


W + jets 


(9516 ± 9) • 10^ 


4.8 ±2.0 


1.6 ± 1.1 


< 1-0 


< 1-0 


di-boson 


(7719 ± 8) • 10^ 


2573 ± 17 


605 ± 11 


56.7 ±4.4 


6.1 ± 1.1 


all SM 


(12540 ± 3) ■ 10^ 


5110 ±49 


1669 ± 32 


503 ± 20 


64.7 ±7.2 


BEl 


23040 ± 47 


14412 ± 37 


13925 ± 37 


12204 ± 34 


11854 ±34 


BE2 


30980 ± 57 


13910 ± 38 


13442 ± 37 


12227 ± 36 


11569 ± 35 


BE3 


31160 ±55 


9118 ±30 


8700 ± 29 


7807 ± 28 


7533 ± 27 



TABLE X: Number of SM background and signal events after each step in the event selection at ^/s = 14 TeV, scaled to 
an integrated luminosity of 10 fb^^. We provide the results for the benchmark scenarios BEl, BE2 and BE3 (Table HI)) . 
The uncertainties correspond to statistical fluctuations. 



(third column of Table |X|, the expeeted SM back- 
ground is reduced to roughly 5110 events. Already at 
this stage, the expected signal event yield of the bench- 
mark points BEl, BE2 and BE3 is overwhelming, i.e. 
a factor of 2-3 larger than the SM backgrounds. The 
signal efficiency of this first cut is the same as for the 
LHC at = 7 TeV. 

The jet multiplicity requirement (fourth column of 
Table |X]) reduces the SM background to 1670 events. 
It mainly originates from Zj (26%), ti (24%) and WZ 
(15%) production. Because sparton pair production 
strongly dominates the signal for all benchmark sce- 
narios at \/s = 14 TeV, cf. Table |Vl almost every 
signal event has at least two hard jets. Therefore, the 
signal efficiency of this cut is large, i.e. 95% or higher 
for all benchmark points. This is higher than for the 
^ = 7 TeV sample, cf. Table Ell 

The Z veto (fifth column of Table |X| effectively re- 
duces the Z + jets and di-boson backgrounds, leaving 
only a total SM background of roughly 500 events. 
The background is now dominated by the tt produc- 



tion. The number of signal events is only reduced by 
roughly 10%. 

Finally, after the requirement on the visible effective 
mass (last column of Table |X]), the SM background 
is reduced to approximately 65 events. At the same 
time, nearly all signal events pass this cut. The signal 
to background ratio is now of 0(100). This justifies 
neglecting the SM background events for the mass 
reconstruction, cf. Sec. |Vl 



Appendix C: Kinematic Endpoints of Invariant 
Mass Distributions 

Assuming the cascade decay of Fig. [131 analytic for- 
mulas for the (measureable) kinematic endpoints of 
the two- and three- particle invariant masses, Eq. (P^ . 
can be derived [l05- ,115i] 



max\2 
T-cb ) 

I'm max \ 2 



=(m2 



2 



'm\){m\ — m\)/m 
2 \/ 2 2 \ / 2 

/ 2 2 \/ 2 2 \ I 2 
(m^ - mc)(mc - mB)/mc, 



^ — 1 ^ 



1) ("^0^£ 



-mgm^)(Ti 
2 5 



y or {niD — rriji) if rrig < niAmjj < and mj^m^ < rrigmi). 
- {ml, - mD^J (ml + m\f{m\ + m\f - \^m\.m%m\ ] / (4m|m^). 



(CI) 
(C2) 
(C3) 

(C4) 



(C5) 



These equations can be solved for the unknown parti- 
cle masses in the decay chain. 
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Appendix D: Three-Body Slepton Decays 



As we have shown in Sec. Ill B 31 some regions of 
the la LSP parameter space exhibit the SUSY mass 
hierarchies 



<M^o 



and 



(Dl) 



(D2) 



where £'j^ is a right-handed non-LSP slepton of the first 
or second generation. In this case, the 3-body decays 




i(k2, A;, 



FIG. 20: Feynman diagram for the 3-body slepton decay 



^ Ft ^ ^ ^ ^ Ft ' 



(D3) 



can be the dominant decay modes of the £'j^ and fi. 
This is for example the case in the benchmark scenario 
BEl (BE2) for the flu (n), cf. Table |VlI] (Table EIIIl). 

In ISAJET7 . 64, that we employ to calculate the 2- 
and 3-body decays of the SUSY particles, the decays in 
Eq. (jD3|) are not implemented, because in most SUSY 
scenarios, the fi is considered to be lighter than the 
other sleptons. 

In this appendix, we fill this gap and calculate the 
missing 3-body slepton decays of Eq. (jD3|) . We show 
the resulting squared matrix elements and give num- 
bers for the respective branching ratios. The phase- 
space integration is performed numerically within 
ISAJET. We use the 2- com ponent spinor techniques 
and notation from Ref. jll6| for the calculation of the 
matrix elements. To our knowledge, the calculation of 
the 3-body decays is not yet given in the literature. 



are denoted by with j — 1, . . . , 4. Using the rules 
and notation of Ref. jll6j , we obtain for the amplitude 



•A^ f ■ *\f ■ \ t i{p-ki)-a 



(D4) 



where Oj = y/2g'Nji^ and the spinor wave functions 

arc yi ~ y(fci,Ai) and x\ = x^{k2,\2)- Squaring the 
amplitude then yields 



\M\^ = Ax\{p-ki) ■ ayiy\{p~ fci) • crx2, (D5) 



with 



\ak\ 



Summing Eq. (|D5P over the spins leads to 



Ai ,A2 



2 2 2; 2 

™13'^23 -P k; 



3J 1 



(D6) 



(D7) 



1. Three-Body Slepton Decay 1'^^ 



We now calculate the 3-body slepton decays — >■ 
Eq. (|D3|) . that are mediated by a virtual neu- 
tralino^'^. Because £r and £^ are sleptons of the first 
two generations, we can neglect contributions propor- 
tional to the (i?-parity conserving) Yukawa couplings. 

The relevant Feynman diagram for the decay — 
l'~t^£~^ is shown in Fig. [20l where the momenta 
(j>,ki,k2,k3) and polarizations (Ai,A2) of the parti- 
cles are indicated. The ncutralino mass eigenstates 



In principle, there are also 3-body decays with virtual 
charginos. However, these decays are negligible due to the 
heavier propagators. Furthermore, the right-handed sleptons 
can not couple to wino-like charginos. 



where 

m23 = (p-fc2)'-(fci+fc3)2, (D8) 

771^3= (p-fci)2 = (fc2 + fc3)'- (D9) 

Here, we have neglected the lepton masses, i.e. 
kl.kl = 0, in Eq. (|D7| . 

We now turn to the decay £'-£'1^. The 

respective Feynman diagram is given in Fig. 1211 The 
amplitude is 



iM = {-ia*){-iaj) 



"ip-k^r^ml„ 



ym , (Dio) 



which leads to the following expression for the total 
amplitude squared: 



\M\' ^By.y^ylyl 



(Dll) 
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The Feynman diagrams for the decay fj^ — > r (.ji 
are given in Fig. [22]and the respective matrix elements 
are^H] 



iMi = {-iaZ){~iaY) 



£*N t i{p-ki)-a 



yi, (DM) 



^ (p-fci)2-TO|o 



44, (D15) 



FIG. 21: Feynman diagram for the 3-body slepton decay 



with 



\aj\ m^o 



(P - fcl)^ - "i|o (p - fcl) 



2 — TO^o 



(D12) 



Summing Eq. (|D10p over the spins, we arrive at 

^ |Xp = B{-ml, - ml,+p' + kl) . (D13) 

Ai , A2 

Here, the proportionality to the neutralino mass, m^o, 
in the amplitude, is due to the helicity flip of the neu- 
tralino in Fig. [5T] 

2. Three-Body Slepton Decay t~ t~ l^t^ 

In this section, we calculate the more complicated 
decays fj~ — >• T~f''=^^. On the one hand, the fi is a 
mixture of the left- and right handed eigenstates. On 
the other hand, we cannot neglect the Yukawa cou- 
plings for the third generation. 



with 



The total amplitude squared is 



(D16) 
(D17) 

(D18) 



3,k=l 



^]al*xl{p- kl) ■ ayiy[{p - ki) ■ uxi 



[a^bl*m^o + al.b'^*m^o]x\{p — ki) ■ ayiXiX2 



+ b'!jbl.*'m,^om^ox\x\xiX2 



(D19) 



where 



J* 



(p-A;i)2-m|o (p-fci)2-m?o 



(D20) 



Summing over the spins of the final state leptons, 
we obtain 



4 

J2 \M\' = [(-^23 - fc?)(-m?3 _ (p2 + fc2 _ ^2^ „ ^2^)(p2 _ ^.2)] 

Ai,A2 j,k=l 

- (aJ6"m^o + albYm^a)mriml^ - kl) + b]bl* m^om^o{p^ + ^1 - m^g - ml^)^, (D21) 



where we have neglected the mass of the first or second are 

generation lepton, i.e. fcn = 0. - , „7„ h,\ . ?^ 

Xj 

We finally calculate the related decay iMu = {iaj ){—iaj ) ^ ^ ^^^^ X2X1. (D23) 

T~ l~ t^, where the the exhibits a helicity flip, 

cf. Fig. [231 The matrix elements for these diagrams The calculation of the squared amplitude is analogous 
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t(^-i,Ai) 



n (p) 




^1 (P) 



^(fca, A, 




FIG. 22: Feynman diagrams for the three-body slepton decay fj^ — >■ r t-^t-n- 
r(fci,Ai) rt(fci,Ai) 



^1 (P) 





FIG. 23: Feynman diagrams for the three-body slepton decay f-^ 



to those for the decay f — r if one changes 



the coefficients alj Vj" 



3. Resulting Branching Ratios 

We now briefly study the new 3-body slepton de- 
cays for the Cfl LSP parameter space in the M1/2 — Mq 
plane. In Fig. [24] we show the same parameter region 
as for the LHC discovery in Fig.[TUl Gray contour lines 
indicate sparticle mass differences (in GeV) that are 
relevant for the three-body slepton decays; see cap- 
tions for more details. 



We show in Fig. 24(a) the branching ratio for the 
decay — > /i~e^ej. The dashed (dotted) gray con- 
tour lines correspond to the mass difference between 
the flu and the LSP (xi)- In the white region, the 
jlfi is heavier than the Xi and decays nearly exclu- 
sively via a 2-body decay into the Xi and a muon. In 
the colored region in Fig. |24(a)[ the flfi is more than 
10 GcV heavier than the Sr LSP. Therefore, there 
is enough phase-space for our decay fl]^ — >■ /i~e~e^ 
at a significant rate. We observe that the branch- 
ing ratio increases with M1/2 and is rather insensi- 
tive to Mq. This increase is due to the competing 
decay /i^ — > ^~t~t^ , Fig. |24(b")| becoming relatively 
less important with increasing Mi/ 2] see the discussion 
below. The decay jl^ — ?• /i~e"'"e^ behaves similarly to 



the decay /i^ — )■ fi~e~e~^, although there are some 
small differences due to the different results for the 
spin-summed squared matrix element, cf. Eq. (|D7p 
and Eq. ((DT3)) . 

The branching ratio of the decay fi^ — > ii^r^f^ 
is shown in Fig. |24(b)[ The decay fi^ ji^r^f^ 
behaves similarly. The dashed (dotted) gray contour 
lines correspond now to the mass difference between 
the fiR and the fi (x?)- For light e^, LSP scenarios, 
i.e. at Ml/2 ~ 380 GeV, the fiR decays with almost 
the same rate into the fi and the cr LSP, because 
both particles are nearly degenerate in mass. How- 
ever, the branching ratio B{jl~j^ — >■ ii~T~fi) decreases 
with increasing M1/2, because the fi mass increases 
more rapidly with M1/2 than the cr mass due to the 
left-handed component of the fi. Therefore, at higher 
values of M1/2 the jiR prefers to decay into the cr 
LSP. 

Wc finally present the branching ratio for the de- 
cay Ti r^e^ejj in Fig. 24(c) The dashed (dotted) 
gray contour lines give the mass difference between 
the fi and the cr LSP (xi)- Since the cr and fi 
are nearly mass degenerate for small M1/2, this decay 
is only kinematically allowed for higher M1/2 values, 
cf. the colored region in Fig. 24(c) Here, the branch- 
ing ratio strongly depends on A/q, i.e. it significantly 
increases with increasing Mq. This is because there is 
also the competing i?-parity violating decay fi — )■ ev^ 
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350 420 490 560 630 700 
Ml/2 [GeV] 



(c) Branching ratio for tlie decay fj~ — r~e~ej. The dotted 
and dashed contour lines correspond to the mass difference 
Mf-^ — M^o and Mf-^ — Mg^^, respectively. 

FIG. 24: Branching ratios for the 3-body slepton decays calculated in Sec. ID II and Sec. ID 21 as a function of M1/2 and 
Mq. The other B3 mSUGRA parameters are Aq = —1250 GeV, tan/? = 5, sgn(^) = + and A23i|gut = 0.045. In the 
white region, the decays are kinematically not allowed or heavily suppressed. 



via A231. Thus, only for scenarios with a low mass decays — > r^e^ej (and — > r^e+e^) become 
difference between the Xi ^'^d fi, i.e. where the Xi important, 
in Fig. [22] and Fig. [23] is nearly on-shell, the 3-body 
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